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A parameter-free and relativistic extension of the RuhrPot meson-baryon model is used to define 
the dominant isoscalar meson-exchange currents. We compute pp-bremsstrahlung observables below 
the it— production threshold using a relativistic hadronic current density that includes impulse, wave 
function re-orthonormalization, meson-recoil, NN creation and annihilation, pn~/ + um^f + prff + 
u)rfi vector-meson decay and NA7(-7r,p) exchange currents. We obtain a good description of the 
available data. The NA7(7r) current is shown to dominate the large two-body contributions and 
closed-form expressions for various non-relativistic approximations are analyzed. An experimental 
sensitivity to the admixture of pseudo-scalar and pseudo-vector admixture of the NN-zr interaction is 
demonstrated. We examine the Lorentz invariance of the NN^NN f-matrices and show a dominantly 
pseudo-vector NN-7T coupling renders impulse approximation calculations without boost operators 
to be essentially exact. Conversely, a similar analysis of the AN ^NN transitions shows that boost 
operators and the two-body NAy wave function re-orthonormalization meson-recoil currents are 
required in NN, AN and AA coupled channel i-matrix applications. The need for additional data 
is stressed. 
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I. INTRODUCTION 

The realization that meson-exchange currents play 
a vital role in the description of the low-energy pp- 
bremsstrahlung observables has consequences which are 
only now coming to be understood. For example, the tra- 
ditional objective of pp-bremsstrahlung investigations, as 
indicated in Fig. 1, centers on the capacity of experiment 
to differentiate the accuracy of the off-shell i-matrices 
that are predicted by a range of model-dependent NN- 
interactions. This is now recognized [0j2| as an exceed- 
ingly difficult task and, at best, is contingent on a com- 
pletely reliable and consistent description of the asso- 
ciated meson-exchange currents. As such, a meaningful 
calculation of the pp-bremsstrahlung observables requires 
knowledge of the meson-baryon form factors, meson- 
exchange currents and the NN-interaction within a fully 
consistent and microscopic effective theory. 

Recognizing the importance of exchange currents in 
pp-bremsstrahlung implies a complete departure form 
the conventional approach to the problem and consider- 
ably changes the nature of such investigation. For many 
years pp-bremsstrahlung was regarded as something of 
a special case in nuclear physics because both meson- 
exchange currents and relativistic effects were expected 
to be small. The principle reason for this expectation 
stems from the fact that gauge invariance demands the 
real photon couples only to conserved currents, so that 
the n-body parts of the complete hadronic current Jr n i 
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FIG. 1. Relationship between the on-shell t- matrix 
in pp-scattering, (p,p) and the off-shell t- matrix in 
pp-bremsstrahlung (p,py). The parameters in the 
NN-interaction V defining t=V+VGt (shown as a bubble) 
are fitted to the (N,N) phase shifts. Bremsstrahlung is usu- 
ally considered in nucleon-pole dominance, where (a) ini- 
tial-, (b) final- and (c) rescattering-interactions are retained, 
but all meson-exchange currents are neglected. Within this 
assumption, bremsstrahlung has been regarded as the best 
means of testing the off-shell f-matrix. However, a sensitiv- 
ity to off-shell effects requires a large photon energy, so that 
G=(E — Ho)" 1 diminishes the dominant nucleon pole contri- 
butions and the exchange currents become important. 

for any given NN-interaction Vnn must satisfy, 



V. J + i [H, J°] 
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To isolate the dominant contributions to the observables, 
it is useful to consider only the static limit, where the 
two-body charge density J™ can be ignored. The isospin 

structure of the exchange currents Jw\ for isovector 
mesons (n,p,..) then reduces to (ti x t 2 ) z , which vanishes 
in isospin conserving processes like pp-bremsstrahlung. 
Relativistic processes can also expected to be small since 
the dominant 7r-exchange contributions to the NN-pair 
creation and annihilation amplitudes share this isospin 
structure. Finally, all NN7 couplings with meson-recoil 
terms can be neglected since they are exactly canceled by 
the corresponding wave function re-orthonormalization 
contributions All of this information suggests that 
a static limit description of pp-bremsstrahlung involves 
only the photon coupling to one of the protons either 
before and/or after (but not during) strong interaction. 
The leading-order exchange currents, according to this 
analysis, begin with the 77 (549 MeV), u (782 MeV) and 
e (975 MeV) iso-scalar mesons, and can therefore be rea- 
sonably neglected. 

However, the above analysis is flawed for several rea- 
sons. Even within the static limit there are purely trans- 
verse currents Jt which automatically satisfy V.J t =0, 
so that current conservation places no constraints on the 
manifestly gauge invariant pwj, uiirj and NA7 exchange 
currents shown in Fig. 2. None of these two-body cur- 
rents can be included simply by introducing the com- 
mutators shown in eq. (1.1), yet they all possess non- 
vanishing isoscalar contributions which can be important 
in pp-bremsstrahlung. In addition, as new experiments 
have succeeded in selecting kinematics that escape the 
consequences of the low-energy theorems Q|| and on- 
shell expansions |6|,0] , they necessarily emphasize dynam- 
ics where the (usually dominant) nucleon-pole contribu- 
tions of Fig. l(a-c) are heavily suppressed by the Greens' 
functions accompanying the highly off-shell i-matrix. As 
such, otherwise less important contributions gain consid- 
erable significance in the observables. This shows that 
the pp-bremsstrahlung dynamics involves much more 
than the off-shell i-matrix and the impulse current, and 
appears to share the complexity of other observables like 
np-bremsstrahlung and n + p ^ d + 7. 

A very long list of pp-bremsstrahlung calculations have 
been reported over the last 45 years. We will make no at- 
tempt to review them all since more recent works || pC| 
already contain appropriate citation and serve to remove 
a number of questionable approximations. A notable ex- 
ception to this trend is found in the very detailed r-space 
calculations reported some 25 years ago by Brown |pl,P2[, 
where the rescattering amplitudes of Fig. 1(c) were re- 
tained and eq. (1.1) was used to constrain the longitu- 
dinal meson-exchange currents. Noteworthy calculations 
since that time have generally been less complete, but 
find their merit in the application of more reliable NN- 
interactions and the exploration of coulomb corrections 
171 and relativistic corrections to the impulse current 

iunra. 
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FIG. 2. Currents included in the present calculations: 
(a) impulse current (b) radiative vector-meson decay currents 
VP7 = pnj + ujtvj + P777 + corpy, (c) NA7 n— and p— ex- 
change currents, (d-e) wave function re-orthonormalization 
and (f) NN7 meson-recoil currents, and (g-i) NN-pair cre- 
ation and annihilation currents. In each of (d-f) and (g-i) we 
show 3 of the 12 time-ordered diagrams with the energy cuts 
represented by dotted lines. None of the exchange currents 
(b-i) can be obtained from eq (1.1). 

Only a few of the most recent pp-bremsstrahlung calcu- 
lations (j],[l^,[l9|,^0| attempt a significant improvement on 
the standard set by Brown, although there are important 
technical differences in these works that we will need to 
consider. 

In the present work we will subject the RuhrPot de- 
scription of meson-baryon interactions to the test of 
reproducing the pp-bremsstrahlung data below the 7T- 
production threshold. The reasons for selecting this ef- 
fective theory are 

• A microscopic definition of the strong form fac- 
tors is available from non-perturbative and self- 
consistent calculation |23|] . The results are com- 
patible with skyrme |24] , non-topological soliton 
|B5j and bag-model |26f| calculations, and moreover, 
with an analysis of experiment p7[ . 

• An NN-interaction model using calculated (not fit- 
ted) form factors has been constructed and gives 
an excellent description of the world data for the 
NN-scattering phase shifts p£l . 

• The extension to include meson-exchange currents 
in the calculation of observables introduces no free 
parameters whatsoever [H,[L5|. 
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• A parameter free-extension of the model to define 
the 3-body interaction has been shown to provide 
an accurate description of the triton bonding en- 
ergy @. 

We indicate the relationships between the form factors, 
the NN interaction and the exchange currents in Fig. 3. 
Although such consistent calculations can (in principle) 
be performed for any effective meson field theory, such 
work has so far only been completed for the RuhrPot de- 
scription and it appears that there are severe difficulties 
in obtaining similar consistency in other models. For ex- 
ample, the RuhrPot form factor calculations have been 
modified to adopt the coupling constants of a conven- 
tional boson-exchange NN interaction and yield results 
p3| which can be accurately parameterized as monopoles 
with typical regularization scales of A ~ 0.8 GeV. While 
consistency demands the use of such form factors, con- 
ventional boson-exchange models require f3(i| artificial 
scales ('cut-offs') of A T > 1.3 GeV and A p ~ 1.8 GeV. 
Such an artificial description of the meson-nucleon ver- 
tices necessarily frustrates any attempt to obtain a real- 
istic description of the meson- exchange currents and the 
3N interaction. 



Electromagnetic 
Form Factors 




Observables 



Local Gauge 
Symmetries 
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FIG. 3. In the RuhrPot effective meson field theory, me- 
son-baryon form factors are calculated non-perturbatively and 
the results are used without adjustment as input for subse- 
quent calculations of electromagnetic form factors, the NN in- 
teraction and the meson-exchange currents. Such consistency 
is necessary, for example, to satisfy gauge invariance (which 
relates the meson-exchange currents to the NN interaction) 
and ensure orthonormality of the wave functions (through in- 
clusion of wave function re-orthonormalization exchange cur- 
rents) . 



Coupled channel i-matrices providing a non- 
perturbative description of all possible transitions be- 
tween NN, AN and AA states have been available 
for more than 25 years l3l|-|33]] and have already been 
used to calculate the A-isobar contributions to pp- 
bremsstrahlung observables 19 2(J. Under these cir- 
cumstances it may appear curious that we choose to 
develop a perturbative description of the NA7 ir— and 
p— exchange currents. However, the coupled channel 
^-matrices used in recent pp-bremsstrahlung calcula- 
tions are obtained by the inconsistent combination of 
the Paris NN^NN and a static limit version of the 



Ried-parameterized Bochum |}2| NN^ AN interaction. 
It is therefore impossible to accurately remove the double 
counted two-pion exchange amplitudes with intermediate 
NA states, so a free parameter is introduced to permit 
an approximate subtraction procedure. As such, this ap- 
proach discards from the outset any hope of a obtaining 
a microscopic description and the quality of the results 
must be interpreted in terms of a meaningless parame- 
ter. We will later show that, even if these inconsistencies 
were to be resolved by fully consistent calculation, such 
an approach is contingent upon a reliable description of 
boost operators, as well as (a subset of) the relativistic 
meson-exchange currents that will be developed in this 
work. 

In the present work we develop our earlier descrip- 
tion 0JTJ,Q of the RuhrPot meson-baryon interactions 
in the pp-bremsstrahlung data below the 7r-production 
threshold. In refs we included the relativistic sin- 
gle and rescattering impulse-current amplitudes, and in 
refs |j],[l5| we introduced the fully relativistic descrip- 
tion of the radiative vector-meson decay currents and the 
non-relativistic description of the NA7 n— and p— meson 
exchange currents without recourse to the soft-photon 
approximation. In the present work we investigate a 
number of important extensions. In particular, after 
describing our model-independent formalism in sec 0, 
we provide in sec III the first bremsstrahlung calcula- 
tions including a fully relativistic description of the wave 
function re-orthonormalization and meson-recoil currents 
that are required to ensure the orthonormality of the 
wave functions is preserved. We also investigate the 
Lorentz structure of the NN-7r vertex by providing the first 
calculations for the purely relativistic NN pair creation 
and annihilation currents. We further present relativistic 
expressions for the NA7 tt— and p— exchange currents 
and identify the source of error in various approxima- 
tions. Supporting calculation details are supplied in the 
three appendices. In sec [V, after establishing the sen- 
sitivity of the selected pp-bremsstrahlung observables to 
each of these currents and concluding that a relativis- 
tic calculation of the isobar amplitudes is necessary, we 
compare our relativistic results with the complete data 
set available from the 1990 TRIUMF pp-bremsstrahlung 
experiment |3f| . We obtain a good description of the ex- 
perimental data and conclude that a large pseudoscalar 
admixture in the NN7r Lagrangian is ruled out. Further 
conclusions and future objectives are given in sec [v|. 



II. FORMALISM 

A. Observables 

We begin by presenting the model independent expres- 
sions we require for the calculation observables for the re- 
action N+N — > N+N+j. The S'-matrix from covariant 
perturbation theory 
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Sf. =S fi -i d A x(f\J^x)A^(x)\i) + 



(2.1) 



gives the probability amplitude for a transition \i) — > \ f) 
as a series involving 0,1,. . . interactions where the pho- 
ton field A^{x) couples to the hadronic current density 
J ll (x). Since only terms with an odd number of electro- 
magnetic interactions can contribute to the production 
of a single real photon, and each of these diminishes by 
a ~ 1/137, we retain only the lowest-order contribution 
and define the transition amplitude as, 

(27r) 4 T fl S^(P f - Pi) = -iJ d A x(f\J^(x)A fl (x)\i) 

(2.2) 



Following a well trodden path, we integrate the transi- 
tion amplitude over the phase space available to the final 
state, and divide by the incident flux, so that with plane 
wave normalized to a (5-function, we obtain the Lab- frame 
differential cross section as, 



d 3 a 



(20 



dfl 3 dfl4d9 1 



M fi | 2 J P 



(2.3) 



where, for the pp- and nn-bremsstrahlung reactions we 
have 



\M% 



M ft \ 2 = l 1 



z E EEi^ 



SiM s . SjM Sf A=l 

and for np-bremsstrahlung we require, 

f \ E i-^i 2 



Mr 12 



Ti,T f 



if , Ms f , A measured 



if , Ms, , A not measured 



if , Ms } , A measured 



a E E E E l^/il 2 if Ms t , Ms f> X not measured 



(2.4) 



(2.5) 



Q 

Ti,Tf SiMs i SfMsj, A=l 



with the invariant amplitude given by, 

M fi = i{2Tr) 15/2 m~ 2 [2uE 1 E 2 E 3 E i ] 1/2 T fl (2.6) 

We retain only the transverse polarization vectors for the 
real photon since, within the Gupta-Blculcr quantiza- 
tion formalism, the longitudinal and scalar components 
can be made to cancel with a gauge transformation, and 
therefore cannot effect the observables. The 'phase-space 



the limit 0~ 



0, where $ — > is guaranteed, N 



factor' J ps appearing in eq (2.3) is defined for arbitrary 
non-coplanarity $ = (ir + 3 — 4>a)/2 as, 



2 2 

j = vivi 
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N 



k sin 9~f cos 9 1 



(2.7) 



with 



N = (p 4 sin 64 - p 3 sin# 3 ) [sin(0 3 + 6a) 

— (f3 3 sin 9 a + (3a sin 9 3 ) cos 7 ] 
— k sin 2 9 1 (/3 3 cos #4 — (3a cos 9 3 ) 
+ 2 sin 9 3 sin 9a sin 2 $ [ (p 3 cos 9 3 — p 4 cos 9a) 

-(p3/3 3 -P4/3 4 ) cos 7 ] (2.8) 

where we use (3i—pi/Ei for laboratory-frame reaction 
kinematics pi + P2 —> P3 + Pa + k. We realize that in 



possesses singularities but ksvc\9 1 /N remains well de- 
fined. Our description of kinematics and phase space 
is the same as that reported in the detailed discussion of 
ref |3(|. As such, it is sufficient to note that in rcf [^6| it 
was shown that J ps possess a square root singularity at 
the kinematic limit of non-coplanarity, although in the 
present work we avoid the non-relativistic simplifications 
discussed therein. 

For the calculation of polarized observables it is con- 
venient to denote da(±i) as the cross section of eq ( |2.3j ) 
measured with the beam polarized in the ±i direction. 
We choose the quantization axis as the beam direction in 
the lab frame, and define the vector analyzing powers as, 



da(+i) - da(-i) _ Er.r.AT^/i^'Oli]^/,} 
da(+l) + da(-i) ET^A^f^/^/i} 



(2.9) 

where i = x,y or z in the lab frame. Similarly, the ten- 
sor analyzing powers (sometimes called 'spin-correlation 
coefficients') are given by, 

_ da(-H,+j) + da(-i-j)-da(-H-j)-da (-4 , +j ) 

t\- % j — ~ ~ ~ ~ ~ ~ ~ ~ — 

da (+i ,+j)+da(-i,-j)+da(-U-j) + d<7(-i, +j ) 
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E 



TiTfX 



(2.10) 



particular, we can expand both A and the Hamiltonian 
H in powers n of the coupling constant, so that with the 
free-particle energy denoted as H , we have, 



where, for example, da(+i, —j) is the cross section mea- 
sured with the beam polarized in the +i direction and 
the target polarized in the — j direction. 



n=l 



and note that rjHirj=riHoX=XHQi]—0, to obtain, 



(2.15) 



B. The Hadronic Current 

To obtain a microscopic definition of the invariant am- 
plitude we compute the Fock-space matrix elements of 
the ph oton field and hadronic current densities appearing 
in eq (2.2), so that after making use of translational in- 
variance and selecting the Lorentz-Heaviside system with 
natural units, eq ( |2.6| ) can be recast as, 

M fi = (2Tr) 6 m- 2 ^E 1 E 2 E 3 E 4 e^k, A)^! J"(0)|^ +) ) 

(2.11) 

where the field-theoretic hadronic current is given by, 



dC 



dC 



d{d v A») OA, 



(2.12) 



for the Lagrangian £ describing electromagnetic inter- 
actions with the interacting meson-baryon system. Di- 
rect calculation of eq. ( p. 11 ) is impossible since |^E'| + ' > ) 



and |^ ■*) represent complete meson-baryon states and 
therefore involve nucleon, resonance and meson de- 
grees of freedom to infinite order. The problem can, 
however, be approached with a Hamilton formalism 
H where the total Hilbert space is partioned into 
meson+resonance+anti-nucleon vacuum- and existing 
subspaces, 

H r , = {iNN)}, Hx = {|the rest)} (2.13) 

We will refer to these subspaces as the ?7-space and A- 
space respectively. Defining r\ and A as operators sat- 
isfying the conventional algebra rj 2 — 77, A 2 = A and 
r/A = Xr] = and which project out the components of H v 
and H\ respectively, we apply a unitary transformation 
to decouple the meson-resonance vacuum and existing 
components of wave functions. Although this formalism 
makes provision for applications involving (for example) 
explicit meson and/or A (see App. A) degrees of free- 
dom in the initial and final states, we confine our present 
application to energies below the 7r-production thresh- 
old, so that the complete interacting meson-baryon wave 
function can be written as, 



1 



\x) 



(2.14) 



where \X) is the two- nucleon state vector and A = XArj 
is required to satisfy X(H + [H, A] - AHA)rj = 0. In 



H n + [H Q , A n ] + ^ HiA n -i 



n=l 

n—2 n — i — 1 



; = 1 



i=l j=l 



1-3 



(2.16) 



We are free to further constrain A by demanding 
eq. ( [2.16] ) is satisfied at each order of n, as would 
be required for any perturbative application. Since 
Ho rj\^/)=£i 7]\^>), where £i is the asymptotic energy of 
the free two-nucleon state, we obtain, 



(Si - H )A n = X 



n-1 

i=l 
n — 2 ri— i — 1 

~ y y .\,n,.\„ , , 

i=l 3=1 



V (2-17) 



Since A o =0, we have 
A 



Ai = 
A 2 = 



£■% — H 
X 

£i — Hq 



H 2 V 



(2.18a) 

A -Hi - X TT H x j] (2.18b) 



£i — Hq £i — Hq 



Finally, with 



Hi 



Cd 3 x 



(2.19) 



we observe that A is completely determined by the strong 
interaction Lagrangian densi ty de fining any model of in- 
terest. Combining eqs. (2.11) and (2.14), we then obtain, 



M fi = {2n) e m^ 2 ^/E 1 E 2 E 3 E 4 (X f \ e u (fc, A) j£ ff (0) \X t ) 



where 

jyo) = ?? ^=L==(i + At)^(o)(i + A) 1 



(2.20) 



n 



effv_/ ''VT+aJa k ' " ' "'VTTaTa 

= v[,P l (0) + J»(0)A + AU^(0) + A+ J»(0)A 



(2.21) 
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is the effective meson- baryon current density. Eq. (2.21) 
provides a time-ordered relativistic description of the 
impulse- and meson-exchange currents implied by the 
strong interaction Lagrangian density defining any model 
of interest. This provides, without approximation, a non- 
covariant three-vector representation where all particles 
are confined to their mass shells and energy need not 
be conserved at individual vertices. An intuitive under- 
standing of the processes embedded in the effective cur- 
rent density can be obtained by noting that the operator 
A is always associated with transitions from the 77— space 
into the A— space, so that with our present definitions, A 
serves to create meson+resonance+anti-nucleon existing 
states and A^ serves to restore purely two-nucleon states. 
In the second-order expansion of this current we observe 
direct terms J, initial- and final-state interaction terms 
J A + A' J, meson- recoil terms A' J A, and wave function 
re-orthonormalization terms JA^ A and A^AJ. The wave 
function re-orthonormalization terms result directly from 
the requi rement that the transformation used to obtain 
eq. (2.14) is unitary - or equivalently, from the fact that 
we insist upon working with orthonormal wave functions. 
This point has been discussed in considerable detail else- 
where 

We select a momentum-space representation and per- 
form a i-matrix expansion of the two-nucleon wave func- 
tions according to the standard procedure. With the pho- 
ton field quantized in the Gupta-Bleuler formalism, we 
require only the transverse polarization vectors. Since 
these have a vanishing time component, we require only 
the spacial parts of the effective current density, so that, 



Mfi = Ne(k, A)(p 3 j?4;a/|Jg ff (0)|pip 2 ;ai) 

(2.22a) 

+ Ne(k, X)(p 3 p^a f \J eS (0)G l t^\p 1 p^a l ) 

(2.22b) 

+ Ne{k, X)(p3P^a f \t^G f J cS (0)\p 1 ]h{ az } 

(2.22c) 

+ Ne(k, X)(p 3 p^a f \t^G f J eS (0)Ga {+ ^PiA-a l )} 

(2.22d) 

where Gi and Gf are ?7-space Green's functions describ- 
ing the propagation of two-nucleon states and 



N = -(2n) 6 m- 2 ^/E 1 E 2 E 3 E i 



(2.23) 



The four terms shown in eq. (2.22) will be referred to 
as 'direct', 'initial-state', 'final-state' and 'rescattering' 
amplitudes respectively. 



III. MODEL DEFINITION AND CALCULATION 
DETAILS 



A. The RuhrPot Lagrangian 



We adopt the strong-interaction Lagrangian densities, 



*-NNtt 

£nN<5 

•£nAtt 
^NAp 

where 



.T 



-«ffNN7rV'[A7 5 - (1 - 

-5nn p '07 a1 V'P/ 
-g-NNSipipS-T 

SNAtt 

2m 

■ ffNAp 



Ip^TNAtpd^TT + ll.C. 



?/^7 5 7 l 'r N A#/ii/ + h.c. 



2 m 



(3.1) 



= d^pv — <9i,Pp, and we have denoted 
.gNAp^MNAffp, ffNN P =5ffp and 5nno; = 55w i where g p and 
g^ the strong charges for the p and w gauge fields. For the 
electromagnetic-interaction we use the Lorentz-Heaviside 
system with natural units, where the charge of the pro- 
ton is defined as e p = +\/4ira with a ^1/137.04, and we 
adopt the Lagrangians, 



£nn 7 = -e N V'7 M '0-4p + -^-ipa^duA^ijj 



£pv 7 
£na 7 



2my 



2m 

-<>'" tT F ia ,$Z.dr< 



where F uv = duA v - d v A 



(1+r 3 ) 



and /lin= 



pi 

(i+« 1B ) 
2 



x. (k 



p or uj, 



(3.2) 



-tt or 77, 



1+kn with 



P 2 f*JN — 2 1 2 

K is =-0.12 and K iv =3.706. 

These Lagrangians describe the bare coupling of mass 
renormalized fields. The form factors describing the cou- 
pling constant renormalization have been calculated |23]j 
as a coupled set of integral equations yielding results 
which can be accurately parameterized as monopoles 
with the regularization scales shown in Table 1. These 
renormalized couplings are implemented by replacing the 
Lagrangians with the vertex functions describing their 
dressed counterparts, 



Tnnu- 

TNNr; 
LNNn = — 



-ig^F NN ^[Xj 5 - (1 - A)-^- 7 5 7 ' i (0 Al )]V'7r.f 

2m 



i -ipj 7 9/<W 

2m 



NNp 
TnN<5 



g^lF^r - 

gNN8F NN 5?p1p5.T 

ffNNc^NNeV'V'e 
SNAtt-FnAit 



KnF 



(2) 



NNp 

2m 



k F (2) 

" NNm „iiv a l„/„ , 

2m 



2m 



^TNAipd^ir + h.c. 



G 



- NAp 



where 



.ffNNpG N Ap ffNATr ju 5 „_ , 

-■0^7 O 7 r NA Vw + h.c. 



2m 

(1) 



(3.3) 



SNNtt 

K p^NAp an< ^ we normalize all form 



r NAp I NAp 

factors as F(0) = 1. As in other exchange current 
applications |37j| the experimentally unknown values of 
5NAp = 20.73 and /iNA=3-993 are fixed according to 
SU(6) p8|,^9| and vector-meson dominance as, 



Mna = ' 



where 



, SNAtt 
5NNtt 

M(0) 



5NAp = 5NNp(l + Kp) 



3NAtt 



(3.4) 



2.353. 



ffNNn 

Note that the tensor cou- 



plingSKp and absent from the Lagrangians of 

eq. (p.l| ) but appear in the dressed vertex functions of 
eq. ( |3.3| ) since they are directly computed from the loop- 
integrals appearing in the form factor calculation p3| . A 
similar consideration for the electromagnetic form factors 
is obviously not required for the real photon. 



The NNa properties are taken from the fit of 
the RuhrPot two-nucleon interaction [Q to the NN- 
scattering data. The form factor scales adopted in 
rcf f2§| l were actually calculated within a non-relativistic 
framework but the recent relativistic calcula- 

tion HjJ has confirmed the parameterizations. We ac- 
knowledge some ambiguity in the signs of the PV7 cou- 
pling constants [|||J4^] but adopt g p7r7 =0.53, g u „ r7 =2.58, 
<7p777=l-39 and <7 WJrr =0.15, as reported in ref [Q. We use 
the experimental result .gNA7r = 28.85, which is consistent 
with the Chew-Low 45 and strong-coupling mod- 
els and, moreover, with the form factor calculations p3| . 
We will not fiddle with these values in order to optimize 
selected experimental results since this would spoil the 
consistency between the calculation of the meson-baryon 
form-factors, the NN-interaction and the exchange cur- 
rents. 





mp (MeV) 


.9NN/3 


Kp (GeV)-* 


S/3 


KV 


SAN/3 


TT 


136.5 


12.922 




49.516 




28.85 


V 


548.8 


6.015 










p 


776.9 


1.651 


6.400 


0.0124 


28.105 


20.73 


LO 


782.4 


4.945 


1.088 


12.379 


0.4334 




s 


983.0 


6.043 










e 


975.0 


10.567 




5.6911 







TABLE I. RuhrPot parameters adopted in the present calculation. All NN-meson form factors are taken from direct 
calculation. For the e meson this requires a meson scale of Ai=0.6 GeV, whereas all other mesons require Ai=0.8 GeV. 
(Further details can be found in [23,28]). We adopt the experimental results Qpwy — 0.53, 2.58, gpryy— 1.39, 3^7)7— 0.15 and 

K 8S — _0.12, n lv =3.706. SU(6) and vector dominance indicate /xna=3.993. 



B. Impulse and Exchange Currents 



We describe here the impu lse an d meson-exchange cur- 
rents Jeff, as required in eq. ( [2.22 ). We adopt a partition 
of Hilbert spaces into meson+resonance+anti-nucleon - 
vacuum and -existing parts, as described in section 0. 
In the present work we confine our attention to leading- 
order exchange currents involving the electromagnetic 
coupling to the NN, NN, PV=p7r, lott, pr\ and lot] and 
NA currents, so that, 



and, denoting the nucleon and A-isobar masses as m and 
toa, we condense our notation with, 



Pik 



E ik = \Jvik + to 2 , 
E^ik = \jvik + m 2 A , 



Pi — k for i=l, 2 
Pi + k for «=3, 4 

£ik = Eik + to, 



£auc = Ea 



ik 



TOA 



(3.7) 



t _ r NN i tNN . tPV i tNA 
•Jed — J C S "+" J c ff ' J cS "+" ^eff 



(3.5) 



The effective cu rren t can th erefo re be derived unambigu- 
ously from eqs. Q, Q, ([111 ), fl2.18H2.19p and (^2l]) . 
Throughout we describe the momenta of a meson with 
mass mp with, 



1. Impulse and Exchange currents with the Relativistic NNj 

vertex 



qi=P3-Pl, Wi: 
q2=P4-P2, W 2 = 



•ml, q° 1 =E 3 -E 1 , Q\ = -q\ 



m% q° 2 =E A -E 2 , Q\ = -q\ 

(3.6) 



Fo r the partition of Hilbert spaces defined in sec- 
all contributions involving a vertex where 



II B 



tion 

the photon couples to the nucleo n cu rrent must satisfy 
AJnn?7 = 77</nnA = 0, so that eq. (|3.5[) requires, 



7 



(a 




d) 




FIG. 4. NN7 wave function re-orthonormalization and 
meson-recoil exchange currents. These currents are necessary 
to preserve the orthonormality of the initial- and final-state 
wave functions. 



D 



abc 



_ 1 

2 

[E3 + E4 — E^k — E-2 — L02] [Ei — E 3 k — W2] 
+ 1 

[E4 — Ei — uj 2 ] [Ei — E 3 k — ^2] 
_ 1 

+ 2 - 

[E3 — Ei — 0J2] [E2 — E4 — UJ2] 



u dcf — 



[E3 — Eik — ^2] [^1 + E2 — Eik — E4 — 012] 
+ 1 

[E 3 — Eik — 1^2] [E 2 — E4 — L02] 



+ 



[E4 — E2 — U>2\ [El — E3 — L1J2] 



In the static limit we note that -D^bc = D h Aei 



(3.10) 



so 



that the wave function renormalization and meson recoil 
exchange currents J^ lrr simply vanish. The same conclu- 
sion can be reached within the soft-photon approxima- 
tion in the barycentric frame. 

Since we will avoid these approximations, we are forced 
to accept that a relativistic description of the pho- 
ton coupling to the positive-frequency components of 
the impulse current Jnn 7 necessarily leads to an effec- 
tive current density J c g comprising both one- and two- 
body operators. These two-body contributions have, 
to date, never been explicitly included in any of the 
bremsstrahlung calculations that seek to include the rel- 
ativistic components of the NN7 vertex. 



2. Pair Currents 



(PsPil [J c T] \P1P2) = Jn^[1]{PuP3)8(P4, - pa) 

+ J NN 7 [ 2 KP2,P4)<S(P3 - Pi) 

+ (fm\J^„\pm) (3.8) 

where the first two terms describe the impulse currents 
for nucleons 1 and 2, and the last term denotes the 
wave function re-orthonormalization and meson recoil 
currents, 



(P3P4\J^[„\P1P2) = -9a 







^ab^ J NN 7 W (P3k , P3)H^ Nf3 [1] (pi , P 3 fc)#NN/3 [ 2 ] (p2,Pi) 

+D^ f P K N0 [1] (pik,Ps)H^ m [2] (p 2 , p4)J^, ( [1] (pi k ,pi) 



+(1,3 #2, 4) 



(3.9) 



where 13 = it, r), p, u>, S and e, and the factor —g ar 
resulting from the contraction of the vector meson polar- 
izations and all references to the Lorentz indices a and r 
are to be ignored for the scalar mesons. Explicit expres- 
sions for the vertex functions HnN/3 an d current Jnn 7 are 
supplied in app. B. The propagator functions are labeled 
in correspondence to Fig. 4 and are defined as, 



The sum of one-body impulse-currents J^n 7 anc ^ the 
wave function re-orthonormalization and meson-recoil 
exchange currents J^ frr do not exhaust the requirements 
needed to obtain a relativistic description of the photon 
coupling to the nucleon current density since the off-shell 
nucleon comprises a linear super-position of positive and 
(so far neglected) negative frequency components. In the 
Feynman-Stuckelberg approach, the negative-frequency 
components of the off-shell nucleon field are interpreted 
as anti-particles, so we are led to introduce the photon 
coupling to the NN-pair creation and annihilation cur- 
rents. 

Within our partition of Hilbert spaces, the photon cou- 
pling to the NN-pair creation and a nnih ilation vertices 
must satisfy r)J^ N r) = 0, so that eq. (3.5) requires, 



(P3P4| 



J, 



NN 



W1P2) = -9or ^[ 



D 



zb/ H NNf) W (Pi ' P3fc)-H"NN/3 I 2 ] (P2 , Pi) ^N 7 M ^ k ' ^ 



+ £ C 3 J NN 7 W &) H Sl*p W (Plfe > &)Eh*(> PI (P2 - P4) 



+(1,3 #2,4) 



(3.11) 
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1 


\ 


/ 


\ 

u 


1 


(h) 





FIG. 5. NN pair creation and annihilation me- 

son-exchange currents. These currents are necessary for a 
relativistic description of the NN7 vertex. 

where f3 = tt, T), p, u), 5 and e, and the factor —g aT and 



all references to the Lorentz indices a and r are to be ig- 
nored for the scalar mesons. Explicit expressions for the 
pair-creation and -annihilation currents Jnn 7 and Jjjn 7 
and all required vertex functions are described in app. B. 

The propagator functions are labeled in correspon- 
dence to Fig. 5 and are defined as, 



D p h = 

abc 



1 



-Ezk — El — UJ 2 ][E4 

1 



E%k — Ei — E2] 



[E4 — E 2 



lo 2 ] [Ei 
1 



E 



3 k 



Ei — E2 



D 



dcf 



[— Esk — Ei — uj 2 ][E 2 — L02 — E 4 ] 
1 

[E 2 — Eik — E 3 — E4] [—Eik — Ei — lo 2 ] 
+ 1 

[E 2 — Eik — E 3 — E4] [E 2 — uj 2 — E 4 ] 
+ 1 

[E4 — lo 2 — E 2 ] [—Eik — E 3 — u 2 ] 



(3.12) 



We will adopt these expressions for our numerical work. 
Nonetheless, it is interesting to consider the correspond- 
ing result under various approximations. If we demand 
energy conservation across the current matrix elements, 
then eq. (3.12) reduces to, 



D 





abc 



-2u 2 



1 



{E 3 + k+E 3k )(qi-ml) 



D 



-2lo 2 



def [Er-k+Eik^qi-ml) 



UJ 2 (2m+k) ui 2 m 

2 spa 1 

u 2 (2m~k) L0 2 m 
(3.13) 



where we also provide the static-limit and soft-photon 
reductions. The corresponding 7r-exchange contribution 
under such approximations are, 



(P3P4| 
{PlPi\ 



7 NN 



7 NN 
'off 



\P1P2) 



W1P2) 



(27r) 6 4m(qf + m\) 

spa — Ae p g^ N7r F^ N7r (q 2 ) 
(27r) 6 4m 2 (gf +m2) 



(l + rf)n.f 2 ri.f 2 (l + rf) 



(2m - k) 



(2m + k) 



2) 



Si{B 2 .q 2 ){iTi xf 2 ) 3 + (1^2) 



(3.14) 



Both of these results scale linearly with the parameter 
A controlling the admixture of ps- and pv-couplings in the 
NN-7T Lagrangian, but the well-known isovector structure 
of the non-relativistic pair currents holds only for soft 
photons. Since the data with which we will compare our 
results was planned to maximize the photon energy, we 
anticipate a non-ne gligib le contribution fom the isoscalar 
components of eq ( [3.11 ). This offers the possibility of 
studying A without the complications of describing the 
many leading-order exchange currents that contribute to 



np-brcmsstrahlung. 



3. PVy Currents 



The oj meson (782.4 MeV) decays asw-> ir°j with an 
8.7% branching ratio and indicates the coupling constant 
of (? W7r7 =2.58. As such, the luttj exchange currents can 
be expected to make a non-trivial contribution to both 
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d 





FIG. 6. VP7 = /97T7 + ljtvj + prpi + Ljrjj exchange cur- 
rents. When energy is conserved across these current matrix 
elements, the time-ordered graphs exactly sum to form the 
corresponding Feynman diagrams. 



the pp- and np-bremsstrahlung observables. Analogous 
arguments indicate that the pir^f contributions will be 
large in np-bremsstrahlung, and perhaps also of some 
lesser importance in pp-bremsstrahlung. 

Our desire to preserve complete consistency between 
the form factors, NN-interaction and the exchange cur- 
rents leads us to introduce all leading-order exchange cur- 
rents describing the photon coupling to the decay of all 
vector-meson mesons present in the form factor and NN- 
interaction calculations. We therefore include the PV7 = 
/97T7 + CJ7T7 + prpf + ujrjj exchange currents as shown in 
Fig. 6. Each of these ver tices s atisfie s rjjp yrj = 0, so tha t 
after making use of eqs. j3^), (pTl2t) and ( |2~2T1 ), eq. (§1|) 
requires, 



(P3Pi\ [Jcl] W1P2) = 



•v/4wp(<fi)wv(<f2) 



(27r) 3 m v [gi — Wp][g| - m 2 



x H NNP [1] (pi , p 3 )e v (q 2 , A v )i?NNV I 2 ] (Pz , P4 ) Jyp 7 (qi , q_i ) 
+(1,3#2,4) (3.15) 

where Jvp 7 is given in App. B. We require only the 3- 
vector current, for which the relativistic form is, 



(P3P4\J7b\PiP2) = - 



epgyp-y Fypj \J (qi Vv (ff 2 ) 
(2n) 3 m v [q 2 - m 2 ][q 2 - m v ] 



x J ff NN p[l](pi,p 3 )|iJSNv[ 2 ](P2,P4)['fi x q 2 ] 

+J?nnv[2](p2,P4) x [(g 2 )o<zi - (91)0^2]} 
+ (1,3 ^ 2,4) 



(3.16) 



We will not resort to the non-relativistic limit, but we 
realize it implies q®=q2=Q : so that we recover the well- 
known result for the emission of a real photon of momen- 
tum fc, 



(P3P4|^f P |PlP2 



«e p gvP73NNV5NNP 



(27r) 6 2m v m-[gf + mp][q% + Tri 



xtfi-qiM x ( f 2 )(f 1 ) P (f 2 )v + (l,3^2,4) (3.17) 

where (ti)p(t2)v = T1.T2, f^, t% and 11 for the p7T7, unrj, 
fyrf) and lorpj currents respectively. 



4- NAy Exchange Currents 

All contributions involving a vertex where the pho- 
ton couples to the NA cur rent must satisfy tjJanV = 
vJnaV = 0, so that eq. (|3.5[) requires, 



(P3Pi\J% A \PiP2) 



ffbc J NA 7 [1] (p 3 fc, Ps)^ AN/3 [1] {pi , P3k)H^ N/3 [2] (p 2 , Pi) 
- D d e { H NA0 [1] (Plfc ) P3) JaNj [1] (Pi , Plfc)#NN/3 1 2 ] (j>2 , Pi) 

-(1,3^2,4) (3.18) 



where /3=7? or p and the factor —g aT resulting from the 
contraction of the p-meson polarizations and all refer- 
ences to the Lorentz indices a and r are to be ignored 
for /3=7r. Explicit express ions for the vertex functions 

are supplied in app. B. 



and currents shown in eq. (3.1 



We ignore all negative frequency resonance contributions. 
The propagator functions are labeled in correspondence 
to Fig. 7 and are defined in analogy to the previous sec- 
tions. Using energy conservation for the current ma- 
trix elements and introducing the A decay width Ta via 
p7|,p8[ E Ak — > E Ak — iT A /2, the propagators reduce to, 



D 



13 _ 2^2 

ai " " {ql - m}){E3 + k — E A3k + iT A /2) 



D - 

^def — (JZ 



2LOO 



{qi - m Z p){E x -k- E Alk + iT A /2) 



(3.19) 



An exact calculation of eq. (3.18) can be achieved with 
the use of the vertex functions and currents given in 
App. B. However, at present there exists some uncer- 
tainty in the coupling constants (?NAp and /iNA > so that 
such a rigorous procedure is of limited interest. We sim- 
plify matters by dropping terms of order p 2 / (E + to) 2 
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(a 



d) 



b) 



(e) 



™1 



/ 



(c) 



V 



(f 



FIG. 7. NA7 7T— and p— meson exchange currents. In 
the RuhrPot model the 7r— exchange contributions represent 
the largest of the two-body currents in pp-bremsstrahlung 
but the p— exchange contributions are very small. When en- 
ergy is conserved across these current matrix elements, the 
time-ordered graphs exactly sum to form the corresponding 
Feynman diagrams. 



beyond leading order - an approximation which does not 
involve any p/m— expansion and should be accurate to 
within a few per-cent at energies below the 7r-production 
threshold. This is surel y ade quate for the first-order Sfi- 



matrix described in sec II A 



We introduce the condensed notation 
V 



El -d - Pik 

~ ! ' Aki — ~ 



Qij = Pik — 



■Pj. 



K-i = Pik 

( Pik 



2EAik — Ei 



C 



rriA 
Pi_ 

\£Aik £i 



Pi, 



(m + tua) 



1,2,3,4. (3.20) 



and proceed to calculate separately the contributions 
from it- and p- exchange, 



?NA/ 



(3.21) 



5. NA-y it— Exchange Currents 

The effective current describing the excitation of inter- 
mediate A-isobar s thr ough 7r-exchange can now be ob- 
tained from eqs. ( |3.18 ) and App. B. For the real photon, 
we require only the spacial current, which takes the form, 

(P3P4\jf s A {^)\piP2) 



Nx5 2 .{V A -V 2 ) 



2r|- 



-(jT!XT 2 ) 3 

x l(/Ci+£i)x Qi 3 + (iaix Qi 3 )x£ 1 
AAax x [(£ x - E{) x Q 13 ] - 2iax[{lCx -Ex). Qis 
■N 3 a 2 . (P4-P2) 



2r|- 



(iTxXT 2 f 

{K. 3 +C 3 )xQ 31 -(iaxxQ 31 )xC 3 
ia x x \{K 3 -E 3 ) x Q 3 x]+2iax[(^3-E 3 ).Q 31 } 



+ (1,3 #2, 4) 
where, 

-iep/XNA3NA7r5NNTr 



(3.22) 



N = 
Nx= 



(27r) 6 72m 2 

jVogAlfcg|ef 

2-EAifc 



1 QEi E 2 E 3 Etx 



2 Fn-NttFj 



NAti 



Utt{c%) 



N 3 = 



N £ A3k D* 



2i?A3fc 



(3.23) 



We will adopt eq. (3.22) for our numerical applications 
and use it to consider the merit of various approximations 
that are conventionally adopted to recover a simplified 
operator structure. 

The first level of approximation involves taking the 
stat ic lim it and ignoring the N-A mass difference in 
eq. (3.20). We call this the vertex static limit approxima- 
tion and note that it is equivalent to casting eq. ( 3.22| ) 
as, 

(p3p4|^ff A (7r)|pip 2 ) 

(Nx + N 3 )^Ml [ 4t 3^ _ x f 2 f{iax x q 2 )\ x k 



- 2(Nx - A^ 3 ) 



2m 
(0*2.92) 



2m 



+ (1,3 # 2,4) 



}(iai x q 2 ) - (ifx x T 2 fq 2 ] x k 

(3.24) 



The substantial simplification results primarily be- 
cause the vertex static limit approximation indicates 
Kx—Cx=— k, and IC 3 ~C 3 =k, so that all operator struc- 
tures involving fCi — Ci immediately vanish. However, 
if we consider the static limit with the more reasonable 
approximation ttia ~ |w, then we find 



K,x +£x ~ --jk- ipi, 



Kx - Cx 

K-z - £3 



3 r 13^ 

--K A pi 

8 24^ 

3 r 13^ 

—k A p 3 

8 24 W 

(3.25) 
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so that the /Cj + Ci terms surviving in the vertex static 
limit are rather well approximated, but the neglected 
fC% — Ci terms are poorly represented. 

The second level of approximation involves casting the 
complete expression in the static limit. This is the ap- 
proximation we adopted in earlier work mH|. At the 
small momentum transfers relevant to the present numer- 
ical application, this complete static limit approximation 
will differ little from the vertex static limit and, if we 
further set Ta = and drop the form factor dependence, 
it is equivalent to casting eq. (3.22) in the simpler form, 



(p3P4\jf{t{ir)\pm) = 



-ie p ^NA3NAir3NNir 

(27r) 6 36m 3 



x(CT 2 .g 2 ) 



J 4(m A - m)q 2 + 2|fc|(ioi x q 2 ) 
\ [(m A -m)a-|*H(gf + m»)' 
(m A ~ m)(i&\ x q 2 ) + 2\k\q 2 



[(mi — m) 2 — \k\ 2 ](q^ + m 2 ) 



+ (1,3 # 2,4) 



(in x f 2 ) | x k 
(3.26) 



The third and final approximation we consider is ob- 
tained by neglecting the \k [-dependence in the complete 
static limit description of the baryon propagators. This 
soft-photon approximation gives the conventional descrip- 
tion |j7j of the NA7(7r) exchange current for the radia- 
tion of a photon of momentum k, 



(P3Pi\J c ff MW1P2) = 



(27r) 6 36m 3 (mA — m) 



4r 3 



2 92 



{ (% + m l) 

+(1,3 #2, 4) 



(in x T 2 f(iSi x q 2 ) 
(q 2 + m 2 ) 



x k 
(3.27) 



We anticipate that this result will differ from the com- 
plete static limit whenever the photon energy is compa- 
rable to the N-A mass difference. 



6. NAy p— Exchange Currents 

Within the soft-photon approximation, it is well known 
p7| that the NA7 p-exchange currents are small com- 
pared to the corresponding 7r-exchange currents, al- 
though the destructive interference between the two 
makes it necessary to include both. In the present work, 
we anticipate the p-exchange contributions to be less im- 
portant than in other works because the RuhrPot model 
suggests a very weak NNp coupling constant. (In par- 
ticular, RuhrPot H uses 3^ Np /47r=0.2169 and k p =6.4, 
whereas (for example) Bonn B Q adopts g^ p / At:=Q.Q2 



and k p =6.1). 

Since the NA7 p-exchange currents are expected to be 
small, we proceed by taking the vertex static limit ap- 
proximation from the outset. In complete analogy to the 
derivation of eq. ( |3.22 ), we obtain, 



jNA 



(p) 



(p;E) 



fNA 



(p;M) 



(3.28) 



where, 



(P 3 P4\J% A (P-Mpip 2 ) = {[4iv4 +) (p 2 + P±) x q 2 + 2iV< Ha x x [{p 2 + p 4 ) x g 2 ]]r 2 3 

[+2N^ ] {p 2 + p 4 ) xq 2 + N^iai x [(p 2 + p A ) x q 2 }] (in x f 2 ) 3 } x k + (1,3 ^ 2,4) 



(3.29) 



and 



(P3P4\J?s A (p; M)\pip 2 ) ={ [4i\^ +) [(ia 2 x q 2 ) x q 2 ] ~ 2N^ > [S x x [(a 2 x q 2 ) x q 2 



r(+) 



C-)l 



+ [2JV& } [(ia 2 x q 2 ) x q 2 ] - N^ 1 [a x x [(a 2 x q 2 ) x q 2 \] ] (in x f 2 ) 3 } x k + (1, 3^2,4) (3.30) 



r(+) 



with, 



r(2) 



N, 



' = m4%(Qt)+K P 4\(Qt)}[D p de f-DU, 



M 



No 



SU(6) e p9NNp G M(°) G M P (Ql 

100m 3 (27r) 6 w p (<f 2 ) 



N^=N FX> p (Ql)[D? ei - 



N^=N fW (Ql)[D» Aei -D- 



P 1 

abcJ ' 



(3.31) 



where we normalize Gj^ Ap (0) = 1 + k p and adopt the propagators D p of eq. (3.19). The c omplet e sta tic limit and 
soft-photon approximations follow in analogy to the procedures used to develop eqs. ( [3.26 ) and ( 3.27 )), the latter 
resulting in, 



&pA^s(p-MpiP2) 



47V (1 + K p ) j 4rf (#t + P2) x q 2 (n x r 2 ) 3 ai x [(p 4 +p 2 ) x q 2 ] 



(toa — to) 



(91 



(% + m V> 



x k 
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(p 3 P4\J% a (p;M)\piP2) 



+ (1-2) 

4AT (1 + k p ) 2 J 4r 2 [(^2 x <f 2 ) x q 2 ] (n x r 2 ) 3 [<?i x [(ict 2 x g 2 ) x 92]] 



(toa — w 
-(1-2) 



(gf + m2) 



(91 



(3.32a) 



x k 



(3.32b) 



where J^J*(p;M) is the conventional result, and 
J^ A (p; E) is an additional piece (resulting from the con- 
vection current part of the NNp vertex) which is usually 
ignored on the basis that it is smaller than (p; M) by 
a factor of about 1+k ~ 7. 



C. Direct, Single- and Rescattering Amplitudes 

A model-independent expression for the complete in- 
variant amplitude was given in eq. fl2.22| ), where we de- 
veloped a decomposition into the four terms describing 
direct, single (i.e. initial- and final-state) and rescat- 
teri ng am plitudes shown in Fig. 1. The left-hand side of 
eq ( 2.22) is d efine d by the model- independent expressions 
of eqs. (2.2) and (2.6), and is Lorentz invariant by con- 
struction. However, since the right-hand side of eq. ( 2. 22] ) 
is determined by a model-dependent calculation, there 
can be no a priori guarantee that each of the direct, sin- 
gle and rescattering amplitudes are individually Lorentz 
invariant unless the wave functions are calculated in a 
manifestly covariant formalism pi] ]. 

Since the wave functions are usually constructed from 
an NN-interaction i-matrix that is defined only in the 
barycentric fra me, and two distinct barycentric frames 
appear in eq. ( 2.22 ), we cast the entire expression into 
the (maximally symmetric) average barycentric frame, 
so that the momenta satisfy, 



A — frame : Pi + P2 — = Q = p% + p>4 + ±k 

(3.33) 

and we acknowledge that a formal solution of the initial-, 
final- and rescattering amplitudes requires the applica- 
tion of boost operators [f50|-[53| . 

In the following we will make u se o f the fact that the 
effective current operator of eq. ( |3.5|) is a totally sym- 
metric under interchange of particles 1 and 2. As such, 
we obtain in an arbitrary frame, 

(P3P4; (sis 2 )S f ; {txh)Tf\J eS \pip 2 ; (sis 2 )5j; [tit 2 )T t ) 
= {PAPr, {s2Si)S f ; {t 2 ti)Tf\J ef f\p 2 pi; (s 2 si)Si; (*2*i)T<) 
= (P4,P3; {sis 2 )S f ; (txt2)Tf\J e {i\p2pi; (sis 2 )Si; (tit 2 )Tj) 

x (_l)(S i +S f +T i+ T / ) (3 34) 

Denoting \a) — \(sis 2 )S; (tit 2 )T) , we define the anti- 
symmeterized states as, 



It is easy to see that parity conservation is consistent 
with the Fermi statistics requirement L + S + T=odd in 
the barycentric frame. 



1. Direct Amplitudes 



The direct amplitudes appearing e q. ([j.22] a) can be 
simplified with use of eqs. ( |3.34 ) and (3.35) to give, 



= Ne(k, \).[(p 3 p4; a f \J cS \pip 2 ; ai) 

-(-l) {S ' +T '\p3P4;a f \J cS \p2Pi;a,)] (3.36) 



whe re N is defined in eq. ( 2.23| ) and J c s is given by 
eq. d3.5| ). However, momentum conservation demands 
that the direct terms cannot involve the 1-body part of 
the effective current density, so in the present numerical 
results we include the photon coupling to 

• NN currents with recoil and wave function renor- 
malization currents. 

• NN creation and annihilation currents 

• p7T7, W7T7, prjj and 01777 exchange currents. 

• AN currents with tt and p exchange. 



as described by eqs. |[|), j3Tl| ), ( |3~16| ), ( p2| ), ( |3~29| ) 
and (|3.30|) 



2. Single- Scattering Amplitudes 

T he i nitial-state interaction amplitudes appearing 
eq. 

(|H b ) 

are given by, 



\P\P2\0i) 



-L {\pxp 2 ;a) - (-l)( s+T >|p 2 p i; a)} (3.35) 



M l fi = N^%,\).(0 i p 4 r,a f \I^{Q)]Git^\^p2\CH) 

= Ne(k, ty-^ff dp' 1 dp' 2 (j)3P4; af\J e &(Q)\tf 1 p' 2 ;a) 

a 

x(^ 2 ;a|G i t(+)|^ r 2 ;a i ) (3.37) 

A formal specification of this amplitude follows by in- 
serting the full effective current density of eq. ( |3.5| ) 
and defining boost the operators needed to cast the t- 
matrix in the initial-state barycentric frame. None of the 
bremsstrahlung calculations known to us has attempted 
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either of these tasks. Instead, the current has been trun- 
cated to include only one-body contributions and boost 
operators are ignored under the assumption that the part 
of the invariant amplitude resulting from initial-state cor- 
relations alone is itself individually Lorentz invariant. 

The first approximation could be removed with a 
straightforward application of the expressions provided 
in earlier sections, but since we anticipate the impulse 
current to be significantly larger than the summed ex- 
change currents, our present numerical applications share 
the conventional approximation of retaining only the im- 
pulse current in the single-scattering amplitudes. The 



second approximation will be justified in section IV A 



where we provide a perturbative analysis that indicates 
the NN7 impulse-current contributions to the single- 
scattering terms are close to invariant under Lorentz 
transformation into the barycentric frames. 

We therefore cast the initial-state correlation ampli- 
tudes into the initial-state barycentric frame, where 



I — frame 



p 1 +p 2 = 0=p 3 +p A + k (3.38) 



so that eq. ( 2.22| b) becomes, with the kinematical nota- 
tion of eq. (|3.7|), 



Mfi = Ne(k, X).^2(a f \J N ^[l](p3k,P3)\a) 

Oi 

(p 3k ;a\t (+) \pi;ai) 
X 2E l -2E 3k + ir, 

+Ne(k, A) .y^(a/l J N n 7 [2] (04 k ,Pi)\a) 

a. 

(p 3 ; a|i (+) \pi; on) 



2Ei - 2E 3 



IT] 



(3.39) 



where we have denoted (0, — 0; af\t^ \p, — 0; Qj) = 
{0 ; a/|^ + ' \0; oti). In complete analogy, we cast the final- 
state correlation amplitudes into the final-state barycen- 
tric frame, where 



F — frame : p\ + 02 — k = = 3 + p\ 

so that eq. ( 2.22| c) becomes, 

( P3;a/|^ ( ~ )t |pifc;a) 

2E 3 - 2E lk + in 



(3.40) 



M% = Ne{k,X).J2- 



x(a\J NN -f[l](pi,pik)\ai) 



+Ne(k,\).J2 



(p 3 ;a f \t ( - K\pi\a) 
2E 3 - 2Ei + ir) 



x(a| J NN7 [2](p2,P2fe)|o:i 



(3.41) 



Simple kinematics establishes that the radiation of a real 
photon implies an off-shell i-matrix, so that we arc free 
to immediately take the limits ij — > in eqs. ( |3.39 ) and 
(PI. 



3. Rescattering Amplitudes 

From the results of impulse approximation calculations 
P, ^0|Jl4| we already know that the impulse contributions 
to the rescattering amplitudes constitute a correction of 
< 15% to the corresponding single-scattering amplitudes, 
so for simplicity we neglect from the outset all 2-body 
currents in the rescattering ampl itudes of eq. (2.22d). 
Hence, in the A-frame of eq. (3.33) we obtain, 



a a' \ 

(0 3 ,0 4 ;af\t ( ~' )f Gf\0 lk ,0 2 - 1 a b )(a'\J NNl [l]{0 1 ,0 lk )\a) 

x (Pi , P2 I a a I Git {+) \p\ , p 2 ; on ) 
-(03, PA] a f |* (_)t G/ \p\ , p 2k ; a b ) (a '|Jnn 7 [2] {p 2 , p 2k ) \a) 



x (p\ , P2 ; a a \G t t (+) \0x ,02\OLi) 



(3.42) 



The initial- and final-state barycentric frames differ by 
the photon momentum, so that no frame can be found 
where both i-matrices are expressed in their barycentric 
frame. We therefore introduce a boost operator \ satis- 

tying mm, 



\0a,0 b )= {l-ix{V)}\+P}\V)+0{l/m i ) 1 

V = {0 a +Pb) (3.43) 



0= 7}{0a ~0b), 



Eq. ( |3.42| ) is manifestly symmetric under interchange of 
particles 1 and 2. However, for comput ation p urpos es, it 
proves convenient to make use of eqs. (3.34), (3.35) and 
Q3.43| ) to formally re-express eq. ( |3.42| ) in terms of the 
photon coupling only to nucleon 1, 



M% = 2N(-1)^+ T ^ Y, tj d0(03 + \k;a f \[l + *x(-^)]i ( - )1 G f [l - i X (~)]W+ {%<*') 

a. a.' 

x(a'\J^[l}(-0-^,-0+^)\a)(+0-^a\[l + i X (+^)}Gi^ (3.44) 

l 
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where, since the i-matrix conserves spin- and isospin, 
the sums over the intermediate-state quantum numbers 
(a=S, M s , T, M T ) are restricted such that S=S i: S'=S f , 
T~Ti, T'—Tf with conserved isospin projection Mt- 
For reasons already indicated in the discussion of the 
single-scattering amplitudes, the boost operators can be 
neglected in the present work. A recipe for perform- 
ing the numerical integration over the pole-structures of 
eq. ( |3.44 ) is given in App. C. 



IV. RESULTS AND DISCUSSION 



A. Impulse Contributions 



In the present work we are primarily interested to 
investigate a consistent calculation of the dominant 
isoscalar meson-exchange currents in pp-bremsstrahhmg. 
An important precursor to this lies in establishing that 
the well-known discrepancy between impulse approxima- 
tion calculations and experimental data cannot be re- 
solved by selecting a different (phase-equivalent) NN- 
interaction. Although a qualitative similarity exists be- 
tween the results of recent impulse approximation cal- 
culations for pp-bremsstrahlung observables |p|- ]T4] , ^6| , p8[ 
the calculation differences are generally not confined to 
the differing NN-potentials. Some of these calculations 
describe the photon coupling to the one-b ody current via 
a Foldy-Wouthuysen transformations [[s| ]R||il||il||T^ or 
direct Pauli reduction [ p~0|J^l] , [T^ ] , whereas others adopt 
a static limit description |JT Some include the rescat- 
tering amplitudes of Fig. lc |)] |ll 14 1^] whereas others 
don't [p|, ^2|jr^Jl^ ] . Further differences are found in the 
use of relativistic or non-relativistic two-nucleon propa- 
gators and/or the application of (guessed) off-shell min- 
imal relativity factors |l(],[ll| ■ 

We avoid all of these uncertainties by presenting the 
results of calculations that are identical apart from the 
definition of the potential used to generate the i-matrix 
elements. We also extend the list of commonly com- 
pared potentials to include the Bonn Q , Paris Q , Ni- 
jmegen and RuhrPot pSfl . For the present compar- 
isons we adopt the impulse approximation, so that the 
effective current density developed in section III B re- 



duces to the sum of (one-body) impulse currents, J e s ~ 

>^NN7[1] + Jnn7[2], as defined in App. B, and is there- 
fore common to all potentials. As such, we retain initial- 
state, final-state and rescattering amplitudes with the 
two-nucleon Green's functions described by the relativis- 
tic Lippmann-Schwinger propagators. Partial-waves are 
summed to J max = 8 and no form of soft-photon approx- 
imation is adopted at any stage. 

In Fig. 8 we compare such calculations with 
the TRIUMF coplanar pp-bremsstrahlung data at 
Slab =280 MeV j35|. The cross section geometries are 




28.0°H 
4 = 12.4° 

I I I I I i i i i 



30 60 90 120 150 1 

7 (deg) 



80 



FIG. 8. Comparison of impulse approximation calcula- 
tions using Ruhr, Nijmegen, Paris and Bonn potentials and 
coplanar pp-bremsstrahlung data at i?i a b=280 MeV. The dif- 
ferences between the model results are smaller than their col- 
lective discrepancy with experiment. 

selected to sample both small and large proton emission 
angles. The analyzing power geometry is selected on 
the basis that it is the result most different from zero, 
and therefore presumably the most reliably measured. 
Some differences exist between our analyzing power re- 
sults and those reported elsewhere | ll]] - primarily due to 
differences in the rescattering calculation, as discussed in 
App. C. The essential conclusion here is that impulse ap- 
proximation calculations using Bonn B, Nijmegen, Paris 
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and RuhrPot wave functions are almost indistinguish- 
able, but exhibit a collective discrepancy with experi- 
ment. Given that the purpose of this experiment was to 
distinguish the predictions of such potentials, the differ- 
ences between theory and experiment are large. 

The final task remaining here is to establish that the 
impulse current contri butio ns to t he si ngle-scattering am- 
plitudes given in eqs (3.39) and (3.41) can be accurately 
described without boost operators. This approximation 
is common to all momentum-space bremsstrahlung cal- 
culations known to us, yet it appears to have never been 
verified. Some authors Jl(J|n],^0| have sought a solu- 
tion to the problem by arguing that the non-relativistic 
t-matrix can be made Lorentz invariant simply by at- 
taching 'minimal relativity' factors ]6l| , so that the 
Lippmann-Schwinger (LS) equation can be cast in a form 
that is apparently identical to the Blankenbecler-Sugar 
(Bbs) equation. However, although both of these integral 
equations describe the NN interaction in ladder approx- 
imation, they are not formally identical because the LS 
kernel is constructed in a time-ordered (non-covariant) 
relativistic framework, whereas the BbS-kernel represents 
one of an infinite number of arbitrary three-dimensional 
reductions of the covariant Bethe-Salpeter (BS) equation. 
As such, the Lippmann-Schwinger i-matrix with minimal 
relativity factors should not be confused with a covariant 
definition of the NN-interaction plfl . The most serious 
flaw in the use of minimal relativity factors is, however, 
that the off-shell factors are completely unknown and 
must be simply guessed p^ . We regard this procedure 
as unreasonably arbitrary, particularly since the guessed 
minimal relativity factors contradict the well known form 
of the two-nucleon boost operators []50] [33| . 

Since the i-matrix is defined in the conventional way as 
t=V+VGt but is available only in the barycentric frame, 
we consider the arbitrary-frame perturbative reduction of 
these amplitudes by constructing a toy-model one-boson 
exchange NN-interaction as, 








I I I I I 



30 60 90 120 150 180 

fl 7 (deg) 

FIG. 9. Initial- and final-state interaction amplitudes in 
impulse approximation when the i-matrix is replaced with a 
toy-model boson exchange potential and evaluated in the A-, 
I-, and F-frames. This justifies the neglect of boost operators 
in the present work. 

need to be reconsidered. 



(Pi 1V1 \V\v\,vi) 

= O ^^#NN/3[l](pl,Pl)A 



1 



— H Q £f — H 



x\H NN p[2}(p 2 ,p 2 )j 1 + (1 ^ 2) 



(4.1) 



where /3=7r,ry,p and u), and H^p[i] is the interaction 
energy for a meson coupling to nucleon ?, as defined 
in App. B. In Fig. 9 we present the results including 
only the initial- and final-state interaction amplitudes 
when the RuhrPot i-matrix is replaced with our toy NN- 
interaction, the complete e xpre ssi on be in g cas t into the 
A-, I- and F-frames of eq. (|||), ( |3~38[ ), pT0| ). We find 
that the neglect of boost operators repre sents an er ror o f 
about 3% or less. We conclude that eq. ( 3.39| ) and (3.41) 
are essentially exact descriptions of the impulse contri- 
butions to the single-scattering amplitudes. As such, the 
application of minimal relativity factors p0[ , pT 20 may 



B. Exchange Currents and the Relativistic NN7 
Vertex 



In section [II B we showed that a relativistic descrip- 



tion of the part of the effective current operator involv- 
ing the photon coupling to the nucleon current com- 
prises not only the impulse current, but also two-body 
contributions from the meson-recoil, wave function re- 
orthonormalization and pair currents. Although rela- 
tivistic corrections to 1-body NN7 currents have received 
considerable attention in recent pp-bremsstrahlung works 
SO "15 1' t ne 2-body contributions remain to be investi- 
gated. 

The wave function re-orthonormalization and meson- 
recoil contributions are no t expected to be individually 
small, but in section III B we recalled 10] that their non- 



relativistic limits cancel exactly, leaving only purely rel- 
ativistic effects in the pp-bremsstrahlung observables. In 



1G 



the present numerical applications we retain these contri- 
butions in the spirit of exploring the relativistic aspects 
of the NN7 currents. 

By far the most interesting aspect of the relativistic 
NN7 currents is found in the different manifestations of 
the NN creation and annihilation pair currents given by 
various Dyson-equivalent chiral Lagrangians. In partic- 
ular, the simplest meson-theoretic Lagrangian satisfying 
PCAC and the chiral commutation relations is the renor- 
malizable cr-model, for which the nucleon-meson interac- 
tions arc of the form, 



C = —gi\)\o + i'j r.7T 



(4.2) 



With a chiral transformation ip — » (1 + 
cxp[i7 5 r.7?')]-0, Weinberg |55|] showed that, for 
a suitably constrained £ and redefined pion fields 
7? = (2m/<?)£, the Lagrangian tranforms to give NN-7T 
and NN7T7T interactions with, 



2m 



9 2 Z 2 
Am 2 



7 5 ^T.df,TT+ J-T. (7? X ^7?) 



(4.3) 



Similarly, the chiral transformation ip — > (1 + 
£ 2 )~^ exp[i(g/2m)r.(7r' x £)]ip has been shown by Wess 
and Zumino J5(| to yield NN-7T and NNtttt interactions 
with, 



C 



ps 



9 

2m 



4m 2 



■ 1 - 9 -2 

IT ' T '.IT 7T 

' 2m 



ip (4.4) 



In an elegent summary of these chiral Lagrangians, Gross 
p7| noted that both £ ps and £ pv give the correct 7r — N 
scattering lengths and that the NNirir interactions have 
p and a quantum numbers. 

Our immediate interest lies in the pseudo-vector (pv) 
and pseudo-scalar (ps) NN-7T couplings, both of which are 
included via the hybrid form of the NN7r Lagrangian we 
adopt in eq. ( |3.3| ) with < A < 1 . It is a trivial exercise to 
show that the 7r-coupling to the positive-frequ ency com- 
ponents of the nucleon current described in eq. ( |3.3| ) is in- 
dependent of A, so that non-relativistic calculations can- 
not differentiate the ps- and pv-couplings. In relativistic 
applications A = is commonly preferred, presumably 
because the presence of the derivative in the pv-coupling 
provides for the gauge coupling of a photon to the NN-7T 
vertex, so that the purely isovector non-relativistic seag- 
ull exchange currents can be included even when the mi- 
croscopic structure of the NN-7T form factor is ignored. 
However, there are no formal arguments to rule out A 7^ 
and the subject is still under investigation J5^,^9|. 

In Fig. 10 we examine the relativistic NN7 cur- 
rents in pp-bremsstrahlung by comparing the results 
of calculations with RuhrPot wave functions which ne- 
glect (IA) or retain (IA+MEXC) the wave function re- 
orthonormalization, meson-recoil and NN pair-creation 
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FIG. 10. RuhrPot results in impulse approximation (IA) 
as in Fig. 8 compared to corresponding results when the wave 
function re-orthonormalization, meson-recoil and NN-pair 
currents (IA+MEXC) are included with pseudo- vector (A=0) 
or pseudoscalar (A=l) NN-7T interactions. A can be determined 
from pp-bremsstrahlung experiments! 

and annihilation currents with purely pv-type (A=0) and 
ps-type (A=l) NN7r couplings. 

The first observation is the most important one: A can 
be determined by pp-bremsstrahlung experiments! This 
is surely the cleanest probe of the Lorentz structure of 
the NN7r vertex yet considered. We reserve our prediction 
for A until we have included the other exchange currents 
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given in section LEI. 

A secondary observation, which we anticipated in 
ref E3, is that the 2-body contributions are small for 
A = 0. In particular, the A=0 NN7 MEXC contribu- 
tions to the cross section at (LEP) 3 =28.O° and (HEP) 
04=27.8° are almost entirely negligible. However, the ef- 
fects are visible in the analyzing powers, as well as in the 
cross section at 03=12.0° and 04=12.4°, the latter being 
reduced at 7 =O° and 180° by 0.14 and 0.17 pb/sr 2 /rad 



(i.e. 5% and 8%) respectively. While these effects are cer- 
tainly small, they are already comparable to the model 
differences shown in Fig. 8. 



C. Radiative Vector-Meson Decay Currents 

En Fig. 11 we compare the results of calculations us- 
ing RuhrPot wave functions and the relativistic impulse 
current which either neglect (IA) or retain (EA+VP7) 
the relativistic pnj, unj, pr/j and W777 exchange cur- 
rents ^,^3|. The RuhrPot contributions are uniformly 
small, although we observe a reduction in the cross sec- 
tion at (LEP) 03=12.0° and (HEP) 4 =12.4° at 7 =2O° 
by about 0.1 pb/sr 2 /rad (i.e. 10%). 

We have used the fully relativistic expressions of 



cq (3.16) for the numerical applications of Fig. 11, but to 
identify the dominant behaviour of these exchange cur- 
rents it is sufficient to consider only the static limit, ne- 
glect the lj- and p-meson mass difference and neglect the 
77- meson contributions altogether. En the non-relativistic 
limit, this allows us to write, 



(P3Pi\J^\piP2) 



9 p^SNNpIJ '(qi, <?2) 
[J{qi,q2) 



■ J(q2,qi)]n.f 2 

^J{q2,qi)]rl 

(4.5) 



where 



^(9*1,92) 



«ep#NN7r(<7l-gi)(<?l x q 2 ) 



(27r) 6 2m v m[gf 



- ml] [q%+m 



(4.6) 



Since ^nntt, dpi^y and are essentially fixed by ex- 

periment, the freedom in the vector-meson decay ex- 
change currents is limited to the model-dependent values 
adopted for the experimentally unknown coupling con- 
stants (?NNp and (?nnw We have adopted the RuhrPot 
NN-interaction values of <7NNp=l-65 and 3nnw=4.95, and 
note that these values are consistent with the broken 
SU(3) requirement / 5nnp = ^ ■ The matrix elements 



of both isopin operators in eq. (4.5) reduce to unity in pp- 
bremsstrahlung, so that the lottj contribution completely 
dominates the vector-meson decay current contributions 
and the corresponding pirj currents are some 15 times 
smaller. 

The magnitude of each PV7 exchange current is, of 
course, dependent on the choice of NN-interaction. For 
example, the Bonn B model requires (?NNp : =3.36 and 
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FIG. 11. RuhrPot results in impulse approximation (IA) 
as in Fig. 8 compared to corresponding results when the 
VP7= piT'y + CJ7T7 + prff + uirj-y exchange currents (IA+VP7) 
are included. The RuhrPot model has relatively weak NNp 
and NNoj couplings, so these exchange currents are small, but 
comparable to the model differences shown in Fig. 8. 

.gNNw = 17.5, so that #nn^/Snn p = 27 ' which severely con- 
tradicts the broken SU(3) prediction of 9. Although the 
p7T7 and lottj exchange currents have never been included 
in Bonn model calculations for bremsstrahlung observ- 
ables, it is easy to see that these currents alone would 
be respectively more that 4 and 12 times larger than the 
corresponding results of the present calculation. 
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FIG. 12. RuhrPot results in impulse approximation (IA) 
as in Fig. 8 compared to corresponding results with relativis- 
tic (rel), vertex-static limit (vsl) and soft-photon approxima- 
tion (spa) NA7 7T— and p— exchange currents. A relativistic 
description of the NA7 exchange currents is necessary. 

D. A-isobar Currents 

In Fig. 12 we compare the results of calculations us- 
ing RuhrPot wave functions and the relativistic impulse 
current which either neglect (IA) or retain (IA+A) the 
relativistic NA7 tt— and exchange currents, as pre- 
scribed by the relativistic [rel], vertex-static limit [vsl] 
and soft-photon approximation [spa] expressions devel- 
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FIG. 13. RuhrPot results in impulse approximation (IA) 
as in Fig. 8 compared to corresponding results with relativis- 
tic NA7 exchange currents including n + p— exchange [n + p] 
and 7T— exchange [n] only. The RuhrPot model has a relatively 
weak NNp coupling, so that the p— exchange contribution to 
the NA7 currents is small. 



oped in section III B . We do not show results for the com- 
plete static limit since these results turn out to be almost 
indistinguishable from the vertex static limit. In Fig. 13 
we decompose the contributions to the relativistic NA7 
currents into 7r— and p— exchange contributions. Collec- 
tively, these figures show that the RuhrPot description 
of the NA7 exchange currents are very large and require 
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FIG. 14. Contributions from individual currents to the 
cross section. These are (IA) = impulse current including ini- 
tial-, final- and rescattering correlations, (NA7) with summed 
7T— and p— exchange, (VP7) = pnj + uivy + prjj + ujrjj, 
(WFRR+Pair) = NN7 wave function re-orthonormalization 
and meson recoil and NN-pair creation and annihilation cur- 
rents with A = 0. (Full) denotes the sum of all these currents 
(with interferences). 



F. Comparison With Experiment 

In Figs. 15-18 we compare results of calculations us- 
ing RuhrPot wave functions and the relativistic impulse 
current which either neglect (IA) or retain (IA+MEXC) 
all of the relativistic meson-exchange currents developed 
in the present work. These include the wave function 
re-orthonormalization and meson-recoil currents, the NN 
creation and annihilation pair-currents (for both (A = 0) 
pv- and (A = 1) ps- NN7r-couplings), the pnj, toit^f, pr/j 
and exchange currents and the NA7 currents with 
7T— and p— exchange contributions. We observe that 
the sensitivity to the admixture of ps (A = 1) and pv 
(A = 0) couplings survives when the currents are consis- 
tently combined. 

The inclusion of the relativistic exchange currents can 
be seen to provide a reasonable description of the cross 
section data in Figs. 15-16, but Figs. 17-18 show that a 
persistent discrepancy with experiment remains for small 
#3 and large #4. The impulse approximation results for 
the cross section data give y^JA&tmsx = 5.8, whereas the 
complete exchange current calculations with ps and pv 
NN-7T couplings giving Xp S /datum = 6.1 and Xpv I 'datum 
= 4.7 respectively. Moreover, Fig. 15-16 shows that 
adopting a ps NN-7T coupling produces structure in the 
cross section that is simply absent in the data. As such, 
the present calculations indicate that the data favours 
A ~ 0, although Figs. 17-18 show that some of the most 
serious discrepancies with experiment cannot be resolved 
in terms of the NN7T coupling alone. 

The importance of our conclusions on the Lorcntz 
structure of the NN-7T vertex are contingent on a reliable 
data set, so we vigorously stress the need for more pre- 
cise measurements of all observables where the ps- and 
pv-couplings give very different results. 



a relativistic description, but that the p— exchange con- 
tribution is comparatively small. 



G. Some problems with non-perturbative 
descriptions 



E. Comparison of Exchange Currents 

In Fig. 14 we examine the individual contributions of 
each of the 1- and 2-body currents developed in the pre- 
vious sections with the complete effective current of the 
present work. We observe that the NA7 exchange cur- 
rents are substantially larger than all other contributions 
and we recall from Fig. 13 that, in the RuhrPot descrip- 
tion, these are completely dominated by the 7r— exchange 
contributions. We observe minor but non-negligible con- 
tributions from the VP7 currents, and we recall from 
sec IV C that the magnitude of these contributions is es- 



sentially determined by the NNo; coupling constant in 
the NN interaction. 



Strong interaction transitions between NN, AN and 
AA states can be described non-perturbatively with a 
coupled channel i-matrix fl3l|-|3l| and indeed have al- 
ready been used to calculate pp-bremsstrahlung observ- 
ables Under these circumstances it may appear 
curious that we have chosen to present a perturbative 
description of such amplitudes as meson-exchange cur- 
rents. There are, however, a number of difficulties in ap- 
plications of these coupled channel i-matrices, the most 
serious of which appears to arise from the inconsistencies 
that exist between the Paris Q NN=^NN and the Ried 
parameterized version of the static-limit Bochum j3l] |3^] 
AN =^NN interaction. In particular, the conflicting defi- 
nitions of the NN7r and NA-7T coupling constants and form 
factors makes it impossible to reliably remove the 
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FIG. 15. Coplanar pp-bremsstrahlung data at i?i a b=280 MeV and 64 = 12.4° compared to RuhrPot calculations including 
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double-counted NN^NN amplitudes involving inter- 
mediate NA states, so that a free parameter is introduced 
to simulate the necessary subtraction. Any specification 
of the two-body currents would suffer a similar ambiguity. 
Finally, the assumed Lorentz invariance of the NN ^ AN 
transition i-matrices remains to be investigated. 

This leads us to consider a generalization of our formal- 
ism to obtain a fully consistent and microscopic descrip- 
tion of the non-perturbative transitions between the NN, 
AN and AA states. More precisely, we will identify two 
minimum requirements of such an approach that appear 
to have been been neglected to date. 



Recalling from sec. II B the freedom to choose any de- 



sired partition of the total Hilbert space, we now mod- 
ify our earlier choice so that the A degrees of freedom 
are included in the 77-space. (Details can be found in 
App A). Within this approach, the leading order contri- 
butions involve not only the NN7 impulse, wave func- 
tion re-orthonormalization and meson-recoil terms, but 
also NA7 initial- and final-state correlations terms and 
NA7 wave function re-orthonormalization and meson- 
recoil terms, as shown in Fig. 19. Our earlier specifi- 
cation of the NN7 one- and two-body currents remains 
unchanged and will not be further discussed here. The 
additional leading-order contributions involving A iso- 
bars are given by, 



+ J N A 7 [1] (P3k , P3)H AN(3 [1] (fl, p 3k )H^ N p [2] (p 2 , Pi) 

+ (1,3^2,4) (4.7) 

where f3—n or p and the factor —g a r and all references 
to the Lorentz indices a and r are to be ignored for the 
scalar mesons. The exact form of D@ is given in app A, 
but here it is sufficient to note the non-rclativistic limit, 



D 



13 



D 







-1 



2^,3(02) {mA —m + up) 

(rriA — m) 2 
Uj3(q2)[mA - m + wpfo)] 



(4.8) 



The first term in the brackets of eq ( |4.8| ) represents 
the NA7 initial- and final-state strong-interaction ampli- 
tudes, and the second term gives the NA7 wave function 
re-orthonormalization and meson-recoil currents. Unlike 
their analogous NN7 terms, where only one-body cur- 
rents survive in the static limit, the NA7 wave func- 
tion re-orthonormalization and meson-recoil currents do 
not vanish in static limit - even for soft photons (cf. 
ref |2C[|). Indeed noting that to a — m ~ 2m w shows that 
these contributions make an important contribution to 
the intermediate-state AN amplitudes in the low energy 
observables. 
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FIG. 19. NA7 wave function re-orthonormalization and 
meson-recoil exchange currents. These currents are necessary 
to preserve the orthonormality of the initial- and final-state 
wave functions described by NN, AN coupled channel transi- 
tion t-matrices. 



Since the coupled channel t-matrices have been calcu- 
lated only in the barycentric frame, it remains to cither 
calculate their boost operators, or demo nstra te that they 
are effectively Lorentz invariant. In sec IV A (see Fig. 9) 
we demonstrated such invariance for a toy-model boson- 
exchange potential and commented on the reliability of 
guessing minimal relativity factors. We anticipate that 
a similar approximate invariance will probably hold for 
the coupled channel NN^NN t-matrices involving inter- 
mediate AN-states. However, in Fig. 20 we observe that 
the corresponding NA7 initial- and final-state interaction 
amp litud es calculated in the average barycentric frame of 
eq. ( 3.33 ) are seriously different from the correspo nding 
resul ts ev aluated in the I- and F-frames of eqs. ( 3. 38] ) 
and ( |3.40] ). This shows that the leading-order contribu- 
tions to the coupled channel AN ^NN and NN^ AN 
i-matrices are poorly approximated under the assump- 
tion that they are Lorentz invariant. 



We conclude that a non-perturbative description of the 
A isobar amplitudes via coupled channel i-matrix calcu- 
lations involves two complications: 
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FIG. 20. Comparison of perturbative NA initial- and 
final-state correlation amplitudes calculated in the (average 
barycentric) A-frame and compared to corresponding results 
that are obtained in the initial- and final-state barycentric 
frames. The discrepancies show the need for boost operators 
and indicate errors of about 20% result when the NN=f^ AN 
i-matrix is assumed to be Lorentz invariant. 



• There are non-vanishing contributions from the 
NA7 wave function re-orthonormalization and 
meson-recoil terms. The neglect of these two-body 
currents contributes serious errors at low energies. 

• The initial- and final-state NA7 interactions are 
poorly approximated in the absence of boost oper- 
ators. The neglect of boost operators contributes 
serious errors at high energies. 

As such, a meaningful specification of non-perturbative 
A contributions requires the calculation of two-body cur- 
rents and boost operators. This has not been recognized 
in the past. Although the first requirement can easily be 
satisfied by retaining a subset of the exchange currents 
we have presented in this work, there are outstanding 
problems that need to be solved if boost operators are to 
be defined beyond (D(l/m 4 ). Such developments stand 



as a challenge for future theoretical work but can only be 
approached within a model providing a consistent and 
microscopic description of all meson-baryon dynamics 



H. Relativistic Effects 

In earlier work jO we presented selected pp- 
bremsstrahlung observables calculated using RuhrPot 
wave functions and the relativistic impulse current. The 
re-scattering contributions of Fig. lc were retained, as 
were the relativistic pit^f, w7T7, pr/j and Lor/j exchange 
currents. We also included the NA7 currents with it— 
and p— exchange in the complete static limit, as defined 
No form of soft-photon approximation was 



in sec III B 5 



adopted at any stage. 

In the present work we have extend the exchange cur- 
rents to include the wave function rc-orthonormalization 
and meson-recoil currents of eq. (JT^) and the NN-pair 
creation and annihilation currents (for both (A = 0) pv- 
(A = 1) ps- NN7r-couplings) of eq. (3.11) 



and 



that 



required for a truly relativistic description of the NN7 
vertex. We have also replaced the complete static limit 
description of the dominant NA7(7r ) ex change current 
with the relativistic up-grade of eq. ( 3.22 ). 

In Fig. 21 we consider the magnitude of these purely 
relativistic effects by comparing with our earlier descrip- 
tions of the meson-exchange currents. We consider here 
only the pv NN7T coupling (A = 0). From Figs 10,12 and 
13 we already know that the largest difference between 
these exchange current results stems from the relativis- 
tic corrections to the NA 7(-7r) co ntrib ution . In addition, 
from Fig. 12 and sections III B 5| and IV D we note that 
this difference results from the neglect of the N— A mass 
difference in the complete static limit exchange current 
operators used in ref. JlJ. Although we find no need to 
change the qualitative conclusions reported in ref [jjj, a 
comparison of pp-bremsstrahlung calculations with ex- 
perimental data near the 7r-production threshold clearly 
requires a relativistic description of the isoscalar meson- 
exchange currents. 



V. CONCLUSIONS 

We have presented a parameter-free and relativistic ex- 
tension of the RuhrPot meson-baryon model to define the 
dominant isoscalar meson-exchange currents. These in- 
clude the first relativistic calculations for the wave func- 
tion re-orthonormalization, meson-recoil and NN-pair 
creation and annihilation currents in pp-bremsstrahlung. 
We also included the fully relativistic pTt~f, unrj, prjj and 
mrj^i currents and a relativistic upgrade to our earlier 
NA7(7r,p) exchange currents p]Jl5[. 

The results of these calculations show that the meson- 
exchange contributions to the pp-bremsstrahlung observ- 
ables below the tt— production threshold are large. As 
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FIG. 21. Selected results from figs 15-18 using the ps NN-zr coupling with (Full) and without (IA) the relativistic exchange 
currents, now compared with the (CSL) results of ref [1]. The CSL results differ from the full ones in that the wave function 
re-orthonormalization, meson recoil and NN-pair currents are neglected entirely and the NA7(-7r) exchange current is taken in 
the complete static limit. The largest numerical difference in the m eson-ex chan ge cu rrent results can be traced to the neglect 
of the N— A mass difference in the CSL operators. See also sections 1IIB5 and [VD for a detailed discussion. 



such, a meaningful interpretation of experiment obvi- 
ously requires a completely consistent description of the 
meson-baryon dynamics defining the NN-interaction, the 
exchange currents and the form factors that they contain. 

Although the wave function re-normalization and 
meson-recoil contributions are well known to cancel in the 
static limit for soft photons, this cancellation is poorly 
satisfied in bremsstrahlung experiments where the kine- 
matics has been contrived to maximize the photon en- 
ergy. As such, retaining these contributions is necessary 
if the orthonormality of the wave functions is to be pre- 



served. Although these two-body currents have never be- 
fore been included in bremsstrahlung calculations, they 
are necessary for a relativistic description of the NN7 
vertex. 

The motivation for developing such a relativistic 
scheme is found in one of the oldest outstanding puz- 
zles in nuclear physics. The Lorentz structure of the 
NN-7T Lagrangian £nntt is universally accepted to com- 
prise an unknown mixture of ps (A = 1) and pv (A = 0) 
couplings. This mixture cannot be distinguished by any 
non-relativistic calculation and, for quite different rea- 
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sons, makes almost no impact on relativistic calcula- 
tions for the Nir scattering lengths. However, within 
our relativistic framework we have shown that the ex- 
isting pp-bremsstrahlung data indicates A is small. This 
is surely the most reliable assessment of the NN7r La- 
grangian available to date. 

We quantified the VP7 exchange current contribu- 
tions and noted their obvious relationship to the NN- 
interaction. While the small vector-meson couplings in 
the RuhrPot NN interaction render these effects no larger 
than the NN-interaction model differences in impulse ap- 
proximation, we noted the necessity to include these cur- 
rents in calculations for models using large NNw and NNp 
couplings (e.g. Bonn B). However, the largest isovec- 
tor exchange current in pp-bremsstrahlung results from 
intermediate-state isobar excitation via 7r-exchange. We 
investigated a series of approximations for this current 
and found that a relativistic description is necessary. 
Given the practical importance of this result, a compact 
closed form expression was provided and the sources of 
error in various approximations were identified. 

We demonstrated our earlier assertion JT^pTsj that the 
NA7 wave function re-orthonormalization and meson- 
recoil contributions do not vanish - even in the static 
limit for soft photons. Their neglect in recent applica- 
tions |l9|,^0| indicates that the the orthonormality of the 
initial and final state wave functions is not preserved. 
In a perturbative analysis we showed that the assumed 
Lorentz invariance of the NN^NN interactions is accu- 
rate to about 3%, but that a similar assumption [ ^9|j20| 
for the NN==^ AN interaction implies unacceptable errors 
of around 20%. We noted that a non-perturbative de- 
velopment of our parameter-free calculations requires in- 
clusion of the NA7 wave function re-orthonormalization 
and meson-recoil exchange currents (for which exact ex- 
pressions were provided) and the application of boost op- 
erators. The need for more precise experimental data has 
been stressed. 



ACKNOWLEDGMENTS 

This work is supported by COSY-KFA Julich 
(41140512), the Deutsche Forschungsgemeinschaft (Ga 
153/11-4) and BMFT. 



APPENDIX A: NON-PERTURBATIVE NA7 
AMPLITUDES 



with projection operators satisfying n = Vi an d 

ViVj = Vi^ij- Denoting an arbitrary operator causing 
transitions from the r^-space to t he r^- space as r\jOr\i = 
Oji, we require (in analogy to eq ( 2.22 )) matrix elements 
of the form, 



[M fi ]n =Ne(k,\) J2 a f \m[l + 4 



x[J ff].y[l + Gjit)jl']rii\pip2;ai 



(A2) 



where we have denoted the effective 1+2-body current 
density as 



[J cS ] 



m J + ja + a^j + aUa 
-Ija^a- Ia^aj + 



(A3) 



First consider the amplitudes involving an effective cur- 
rent describing NN— ^NN transitions, which must there- 
fore be taken between inital- and final-state wave func- 
tions constructed from the NN— +NN i-matrix. These am- 
plitudes include the dominant NN7 impulse currents, as 
well as contributions from the NN7 wave function re- 
orthonormalization and meson-recoil currents. (Adopt- 
ing (Al) does not alt er t he ex act e xpressions for these 
currents given by eq. ( |3.9| ) and ( 3.10] )). However, we now 
have additional contributions, as shown in Fig. 19, which 
are given by eq. Q4.7D with 



1 



[£3 - Eaxh - up(q 2 )][E 2 - Ei - ujp(q 2 )] 
1 

~2[E 4 -E 2 - M92)][Ei -E 3 - Ljp(q 2 )} 
1 



2[E 3 - E A i k - wp{q%)\ 
1 



D 







[Ex + E 2 - E Alk -E A - 0Jp(q 2 )] 
1 

[E4 - E 2 - LO f j(q 2 )][Ei - E A3k - up(q 2 )] 
1 

" 2[E 3 -Ex- ujf3(q 2 )][E 2 - E 4 - u; p {q 2 )} 
1 

2 [E 3 + E 4 - E A3k -E 2 - u p (q 2 )] 
1 



[Ex - E A3k - 073(92)] 



(A4) 



Within the formalism of sec II B , a non-perturbative 



description of the A isobar contributions requires our 
partition of the Hilbert space to modified such that, 



T~lr)2 ~ 



{|AN>}, 

W n « = {|AA)} > H x = {|therest)} 



H m = {|NN)}, 



H m = {|NA)}, 



(Al) 



where P—tt or p. Unlike their analogous NN7 terms, 
these NA7 wave function re-orthonormalization and 
meson-recoil currents do not vanish in the static limit. 
Their inclusion is required to preserve the orthonormal- 
ity of the hadronic wave functions. 

Next consider the amplitudes involving an effective 
current describing AN ^NN transitions, which must 
therefore be taken between inital- and final-state wave 
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functions constructed from AN 
matrices. 



[M 



-NN and NN^NN t- 



= Ne{k, A).(p 3 p 4 ;a/|»?i[l + 4i )t <?ii] 
t (-)t/ 



X[J e ff]l2[G 2 l*2l' 

7Ve(fc, X).(p 3 p 4 ; a f \r]i[t\ 2 n Gi2\ 



+ (1,3^2,3) 



(A5) 



These are the NA7 initial- and final-state interaction am- 
plitudes and should obviously be specified in a consistent 
frame. However, the i-matrices are available only in the 
barycentric frame, yet are re quir ed for the initial- and 
final-state interactions of eq. (A5) in barycentric frames 
which differ by the photon momentum. The procedure 



adopted in recent works |jlj,|2(| is to attach (guessed) off- 
shell minimal relativity factors and assume this renders 
the initial- and final-state interaction terms individually 
Lorentz invariant. Rather than adopt this assumption, 
we provide an exact specification of the lea ding-order 
contributions to these amplitudes via eq. (4/7) with, 



D 







2(E 3 + E4 — E/\ik — E 2 ) 
1 



1 



(E 3 - EAik - up) (Ei -E% — top) 

> 



2 (Si + E% — £?A3fc — E4, 

I 



(Ei — i?A3fe — Up) (E2 — E4 — LOp ' 



(A6) 



APPENDIX B: VERTEX FUNCTIONS 
1. Interaction energies for strong vertices 



For a meson with momentum q and mass mp we define, 



Nfj = 3NN/3 



{2Tr) 3 8ujpE f E l 



1/2 



(Bl) 



The interaction energies (0|6(p/) — J d 3 xT^p b'(pi)\0) for the coupling of mesons to the positive frequency nucleon 
current are, 



-^NNtt — Aw^NNtt T, i?NNp — ^P^^NNp ^' #NN<5 = NgflNNS T, 

EL _ EL 
£f £i. 



iNNi 



h° 

"NNp 



Ft 



(i) 



K D F, 



(2) 



NNp ~ n 'P ± NNp 



1 Ef+Ej 
2m 



+ 



F (1) +k F {2) 

NNp ^ ft P r NNp 



1 



Ef+Ei 



2m 



Pf-Pi , io-(pf x pi) 



ZlNNp 
/lNN5 = ^NN<5 



+ k F {2) 

NNp t «-p^NNp 



t f ti 



'Pf 






[ £f 







k F {2) 

P NNp ,-, _ n. 



j _ Pf-Pi _ iv-{Pf x k) 
£f£i £f£i 



1 _ g/ji zg.(p) x pi) 



The interaction energies (0|6(py) — J d 3 xT^p d(—pi)\0) for meson couplings to the pair-creation current are, 

N^h^^T, -HfjNp — "'Y^P^NNp ^' ^NN(5 = Nshfj NS T 



-^NNir 



^NNtt 



— +ii*NNir < 



l_(l_ A )^ll + h+(i-A)^ 
m J L m 



NNp 



l NNp 



S (1) +n F {2) 

NNp 1^ P NNp 



NNp 1" fl P r NNp 



2m 



E, 



m 



v-Pf 

£f 

r(2) 



NNp t «-p^NNp 



P/ X p % 

£ f £i j 

"1 % ~ E i 



2m 



a -Pi 



_ K P F NN P 

a+ ^2m- {Pf+Pt) 



o-Pf °jPi 
£f Ei 



F 



(i) 



K Ft 



(2) 



NNp ' "-P 1 NNp 



Ei 

1 + — 

TO 



i(PfXpi) . Sipf.pi) Pf(a.pi) (a.pf)pi 



EfEi 



EfEi 



£ f £ r 



£ f Ei 



(B2a) 
(B2b) 

(B2c) 

(B2d) 
(B2e) 

(B3a) 
(B3b) 

(B3c) 
(B3d) 
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-Ft 



NNt5 



Si 



£ 



f 



(B3e) 



The interaction energies (0|eP(— pf) — J cPTnn/s w(pi)\Q) for meson couplings to the pair-annihilation current are 



'NNtt 



= -iE 



NNti 



l-(l-A)^-! + [l + (l-A)^- 

TO J L TO . 



Pf-Pi 



EfEi 



h° - 

NNp 



'NNp 



NNp ' 


- k F {2) 

P NNp 


1 - 


F (1) H 

NNp " 


h 'P r NNp 


1 - 


NNp ' 




1 + 


-Fnns 


'a.pf a. pi 





E t - Ei 



2to 



TO 

El 

TO 



(i) 



NNp ' 



KnF, 



(2) 



Pf X p t 



(2) 



NNp / -, 



2to 



(p) + Pi) 



P- 1 NNp x ' 



2to 



£ 



f 



£; 



i{Pf*Pi) . viPf-Pi) Pf(v-Pi) (v-Pf)Pi 



£ f £i 



EfSi 



SfSi 



£ f £i 



where all strong form factors are evaluated at Q 2 
For the NA/3 vertices we define ( see eq. (3.7)), 



N, 



A _ WNAttFnAtt 

2m 



£a£ 



{2tt) 3 8lo^E a E 



5NNp(-fM 9 AN 77 



2m 



<7NNti 



£a£ 



(2tt) 3 8uj p E a E 



The interaction energies (0|6a(P) ~ / cPxTan b' (p) | Q) for meson couplings to the N — > A current are, 

(A; 



Hantt — N v /ianttTNAj Han p — N^h ANp e u TNA 



h-ANn 
^ANp 



-cna- 



PA 



P 



g (PA.P) \ 

V«t.a £aitia J 

-PA-&NA . -, , (ctnA-Pa)"" 



S CNA + PA" „ 
TOA CA"1A 

x [(pa-p)p^-(pa-p)a^)] 



(T.pA, , ,CT.p\ _ (v-PA)v(v-P) 

( ^a" ) + ( — + £A£ 



The interaction energies (0|6(p) — f d 3 2;rAN/3&A(PA)|0) for meson couplings to the A — > N current are, 

HnAtt — A^/lNA7r%AJ -£^NAp = ^' t NA/i''NA 



'NAtt — 



"NAp 







^NA- 


PA ^ 



£7 PA-P) 



TOA £a™A 



ff.pA, . t S.p _ (a. p A ) a {a. p) 
{ ^a~ ) + { — );<T + 5a5 



-g N A-PA ?t , rf (^NA-PA) 
i^NA +PA 



TOA 



£a"ja 



x [(p-Pa)p^-(p-Pa)a^)] 



(B4a) 
(B4b) 

(B4c) 

(B4d) 
(B4e) 



(B5) 

(B6a) 
(B6b) 

(B6c) 

(B7a) 
(B7b) 

(B7c) 



2. Currents for electromagnetic vertices 



We define, 



^7 = 7^ { \ (B8) 



1/2 



7 (2tt) 3 {iEfEi 

The photon coupling to the positive frequency nucleon current (0|6(p/) J e g(0)b' (pi)\0) is given by, 
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JNN7 



JNN7 



7 



F 



(i) 



nn 7 
(i) 



kR 



(2) 



NN7 
(2) " 



E f + F 



NN7 ' ftr NN 7 



2m 



tf Ci 



NN7 ' NN7 



P/ _ Pi 



kF, 



Ef + E, 



NN7 



2m 



2m 



{Pf+Pi) 



Pf-Pi , iv-(Pf x Pi) 
p/.pi ia.(p f xpi) 



S f Si 



£f£i 



The photon coupling to the pair-creation current (0\b(pf) J e gt(0)d(— pi)|0) is given by, 



^NN7 



JNN7 



NN7 + ft - r NN7 



NN7 T f "- r NN7 



r NN 7 "T" Kr NNj 



Ef-Ei 



2m 



£ f 



R 



(i) 



NN 7 



Ei 
rn 



Ei 

1 + — 

m . 



«F, 



(2) 



NN7 



2m 



(Pf + Pi 



(- kF (2) 

t «-^NN7 


1 - 


'a.pf 




I £f 


Ei 



Et — Ei 



2m 



o-Pi 

£i 



i(ff x Pi) . v(jPf-Pi) Pf(#-Pi) (v-Pf)Pi 



£f£t £}£% £f£i £f£i 

The photon coupling to the pair-annihilation current (0|d' (— pf) J e ff (0)&^(pi)|0) is given by, 

JnN 7 = -^7 £ A'JnN7 



Jnn 7 



JNN7 



F-r 



(1) 



NN7 



kF 



(2) 



NN7 



NN7 T f "- r NN7 



F, 



NN7 



kF 



(2) 



NN7 



F 7 -F 



2m 



m 



cr-Pi _ 
£* 
kF (2) 

ft " r NN7 

2m 



F 



(1) 



NN7 



kF, 



(2) 



NN7 



1 



E f - Ei 



2m 



°-P! 



{Pf+Pi) 



£f £1 



i(Pf x Pi) . &(Pf-Pi) Pf{v-Pi) (v-Pf)Pi 



£f£i £f£i 

where the electromagnetic form form factors are defined as 



NN7 ~~ 2 NN 7 2 NN7 



1 



-,(l);is 



1 



-,(l);iu 



., p(2) * ,,is pv 

ft N r NN7 _ 2 NN 7 



£f£i 



.(2); 



iti z:i(2);iu 



— K F 



NN7 



=>■ protons 
neutrons 



4n 7 (^ = 0) + « N F^ 7 (fc^ = 0) = (1) + ( 
^ 7 (* 8 = 0) + «NF^ 7 (fc 2 = 0) = (0) + ( 



(2) 



. K 



K" 



2 

K lS — K V 



For the NA7 vertices we define ( see eq. (3.7)) 
N A 



3p G™i <7AN7T 



£a£ 



4FaF 



GJJ(0) = 1 + = 4.706 



7 (27r) 3 2m 2 5nn , 
The photon coupling to the N — > A current (0|&A(P)^cff (0)& t (p)|0) is given by. 
Jan-/ = N A j^ Nl e„(T NA ) 

(ctna-Pa)""' 



7 JAIN7 

ctna-Pa 



; CNA +pa- 



,a.p A , . ,o.p _ 
( ^a" ) + ( — );a 



W1A " ^AWlA 

x [(pa-p)^-(pa-p)a^)] 
The photon coupling to the A — > N current (0|6(p) J e ff(0)& A (pA)|0) is given by, 



{o.pA)S{5.p) 
£a£ 



•/nA 7 — JNA 7 e i'( r NA) 



JNAi 



(7.pA A . ^ (3.pA)a{o.p) 
{ ^A~ ) + { — );a+ £A£ 



x [(P-PAV^-(P-PA)A^)] 



-PA-C NA t 



m A 



_ , (4a-Pa) 
'^ +Pa ^a" 



(B9a) 
(B9b) 

(B9c) 

(BlOa) 
(BlOb) 

(BlOc) 

(Blla) 
(Bllb) 



(Bile) 

(B12a) 
(B12b) 
(B12c) 

(B13) 

(B14a) 

(B14b) 
(B15a) 

(B15b) 
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The photon coupling to the vector-meson-decay current is given by, 

4W><"-."> = ^"Zf^ (2^^ £ " { " £i(fv - Av)( " p) - (w)t (B16) 

where i*yp 7 (fc 2 = 0)=1, qy and qp are the 4-momenta delivered to the nucleons by the vector and pseudo-scalar 
mesons respectively. 

APPENDIX C: THE MOMENTUM-SPACE NN7 RESCATTERING CALCULATION 

The NN7 rescattering amplitudes in bremsstrahlung were first calculated by Brown M] in r-spacc. For p-space 



calculations, a recipe has already been given in ref |10|, although the residue terms are valid only for soft photons and 
the dominant final-state interaction is scaled by (guessed) off-shell minimal relativity factors and cast into the I-frame 
with a barycentric momentum k. This choice of frame does not minimize the effect of the neglected boost operators 
and the inclusion of minimal relativity factors disturbs the convergen ce pr operties of the p-space integral 



Casting the entire rescattering amplitude into the A-frame of eq ( 3.33| ) requires a somewhat different numerical 



procedure. We start from eq. (3.44) and simplify our notation by defining, 



M% = Kff H(p)d6d(f> Hip) = J cJg^P (CD 



with, 



F(p)=p 2 S m6[E(p 3 + lk) + E(p + ±k)][E(ft. - + E{p - H)] 

x J2 & + i*? S f M Sf ;T f M T \[l + i x (_^)]t(-)t[i - i x (~)]\p + ~fc; S f M' s ; T f M T ) 

M S M' S 

x(S f M / s ;T f M T \J NNl [l](-p + ^SiMg-TiMT) 

x(+p- h-S l M s -,T % M T \[l + ix{+^)]G^ + \\-t X (+^)]\p 1 - ^SiMsjTiMr) 
Gi(p)= {pi - \kf -{p- Jfc) 2 + in = A? - (p - cosi?) 2 + in 
G f (p)= (p 3 + \kf - (p + \k) 2 + i Vf = A}-(p+ -\k\ cos i?) 2 + i Vf 

A- = Q|fc|cOS^ — pi .p2 A 2 = Q|fc|cOS^ 

K= ^(-i)(S/+Si+r /+ ro (C2) 

with cos 1? — p.k = sin 9 p sin 9k cos(0 p — 4>k) + cos 9 p cos 9k ■ Within the (maximally symmetric) A-frame the effects of 
the boost operators xi+^k) and x{~\k) can be expected to substantially cancel and will therefore be neglected in 
the present numerical applications. Their inclusion would change nothing that we will discuss in this appendix. 

For A^ A/ < 0, there are no poles on the real axis, but it is easy to see that there are singularies in the Greens' 
functions at p = pi, p%i P3> Pa as wen as (Pi ~ P)-{P2 — p) and (ps — P)-(pi ~ P), the latter two requiring particular 
attention when cosi?i = ±4v / p*Tp2/|fc| or cosi?/ = ±4\/p 3 .p4/\k\ cause Ai — or Af — 0, so that G^ 1 and/or GJ 1 
each have a 2 nd -order pole (on the real axis). We adopt the utilatarian attitude of noting that for experiments below 
the 7r-production threshold, the second-order poles only occur at 9\ ah ,9\ ah < 6°, and even then, only for #^ ab < 1°. 
Since no data exists in this region, we simply defer a treatment of second order singularities and confine our attention 
to the kinematics containing the simple and separable poles, 

p=p ( i ±) = + ^\k\costf± A, if A 2 > then A; = |p 4 (±) - i|fc| costf| 

p=pf z) =-i|£|cos0± A f if A^>0 then A f = \pf } + ~ \k\ cos d\ (C3) 
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and evaluate A? and A 2 to determine if poles exist on the positive real p— axis. Defining such poles to be vectors so 
that pf=p^=p, we obtain, 



H{p) = V.P. / dp 



F(p) 



- Z 



where 



Z = 



iirF(p 



A? - (p - j\k\ costf) 2 A 2 - (p + £|fc| costf) 2 



(C4) 



(+) 



2A i G / (p< +) ) 2A i G / ( P rO 2A f G i (jpf) 2A f G i (pf>) 



,(-) 



(C5) 



remains well defined. This formally completes the specification of the integral. However, for practical purposes, we 
need to add special forms of zero to smooth the divergences near the poles. With a simple generalization of 



V.P. 



1 



1 j ,a 2 + k k — oi . 



r dp 



(C6) 



we obtain, 
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,62 x 
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rfl 
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dp + 



1 



Gf(P) G/(pH)J 2A iG/ (pH) [ (62-p< +) ) (61 -p. 



In 



(&2 -p,H) (61- p|+ 
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(+ 
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In 



(+ 



(d2-p^ +) ) (di — pjr 



17T 



ITT 



(C8b) 



(C8c) 



(C8d) 



where the domains [al:a2] + [bl:b2] + [cl:c2] + [dl:d2] span [0:oo] and contain the poles p\ + \ p\ \ p^ and p^ 
respectively. In the event that p^ does not exist in [0:oo], then our expressions require F(p^) = 0. 
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